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Abstract 

This paper establishes dynamical localization properties of certain families of unitary 
random operators on the d-dimensional lattice in various regimes. These operators are 
generalizations of one-dimensional physical models of quantum transport and draw their 
name from the analogy with the discrete Anderson model of solid state physics. They 
consist in a product of a deterministic unitary operator and a random unitary operator. 
The deterministic operator has a band structure, is absolutely continuous and plays the role 
of the discrete Laplacian. The random operator is diagonal with elements given by i.i.d. 
random phases distributed according to some absolutely continuous measure and plays the 
role of the random potential. In dimension one, these operators belong to the family of 
CMV-matrices in the theory of orthogonal polynomials on the unit circle. 

We implement the method of Aizenman-Molchanov to prove exponential decay of the 
fractional moments of the Green function for the unitary Anderson model in the following 
three regimes: In any dimension, throughout the spectrum at large disorder and near 
the band edges at arbitrary disorder and, in dimension one, throughout the spectrum at 
arbitrary disorder. We also prove that exponential decay of fractional moments of the 
Green function implies dynamical localization, which in turn implies spectral localization. 

These results complete the analogy with the self- adjoint case where dynamical localiza- 
tion is known to be true in the same three regimes. 
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1 Introduction 



The spectral theory of Schrodinger operators and other selfadjoint operators H used to 
model hamiltonians of quantum mechanical systems has a long history. It can be argued 
that on physical grounds the main motivation for studying spectral properties is their close 
connection (e.g. via the RAGE-Theorem) with dynamical properties of the correspond- 
ing time evolution e~*^*, i.e. the propagation of wave packets under the time-dependent 
Schrodinger equation iip^t) = H'ip{t). 

However, the dynamical information following from spectral properties is not very ac- 
curate and examples have been found where spectral properties are misleading about the 
dynamics. In particular, this is the case for operators with singular continuous spectrum 
or dense point spectrum, spectral types quite common in quantum mechanical models of 
disordered media such as quasiperiodic or random Schrodinger operators, see for instance 
[20], p], [29]. 

As a consequence, much of the recent work on hamiltonians governing disordered sys- 
tems has focused on directly establishing dynamical properties. For example, it has been 
shown that Anderson-type random hamiltonians exhibit dynamical localization in various 
energy regimes, the property that wave packets tp{t) stay localized in space for all times, 
see [1], [2], [31], [23], [M], for example. 

The central object of interest in understanding dynamics is the unitary group U{t) = 
^-iHt^ rather than the hamiltonian H itself. As short time fluctuations will generally not 
have a major impact on long time dynamics, one may discretize time by choosing a time 
unit T > and study U{nT) = = as n — > oo, with the fixed unitary propagator 

U = U{T). 

To further stress the role of the propagator as the central object in studies of dynamics, 
consider hamiltonians H{t) which depend periodically on time, H(t + T) = H(t) for some 
T > and all t. In this case the large time behavior of solutions of iip'{t) = H(t)'ip{t) 
is governed by the unitary propagator U{nT,0) = ?7", where U = U{T,0) is commonly 
referred to as the monodromy operator of the time-periodic system. Note that in this 
case U does not have a meaningful representation of the form e*"^ any more. While such 
representations with selfadjoint operators A exist for abstract reasons, the operator A may 
have little to do with the time-dependent hamiltonian H{t). Actually, in the periodic 
or quasi-periodic cases, the operator A is linked to the so called quasienergy or Floquet 
operator [37], [57], [TO], [39]. 

One consequence of this last fact is that it becomes legitimate to study time periodic 
systems by directly modeling the monodromy operator U based on physical properties. For 
example, time dependent analyses of electronic transport in disordered metallic rings have 
been considered within such a framework, [44], [6], [8], [7], [TT]. Kicked systems, often used 
in the study of quantum chaos, provide another example of such models, see e.g. [18], [2T] . 
|22j . |24j . |26j . [12j . [45], |48j . Similar studies were performed on the dynamical properties 
of pulsed systems, given by smooth Floquet operators, in [9], [38], [39], [28] . 

Our central goal here is to investigate the dynamics of one such model , which we call 
the unitary Anderson model, indicating the presence of disorder by the random parameter 
u}. The name is chosen by analogy to the selfadjoint Anderson model, which, in its discrete 
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version on ^^(Z'^), takes the form 

K = ho + V^. (1.1) 

The potential in (jl.ip is given by real-valued i.i.d. random variables {VLu{k)}j^^^d and /iq 
is a deterministic selfadjoint operator, most commonly the discrete Laplacian on Z"'. By 
comparison, following our previous works [40[ I4H [35], for the unitary Anderson model we 
choose a unitary operator on ^^(Z'^) of the form 

= D^S. (1.2) 

Here 5" is a deterministic unitary operator and a multiplication operator by random 
phases, i.e. for every (p G f^(Z'^) and k € Z"^, 

{D^<i)){k)=e-''^^m, (1-3) 

with i.i.d. random phases 0^ taking values in T := M/27rZ. Note that, despite the formal 
analogy, there is no simple relation between selfadjoint and unitary Anderson models. In 
particular, due to non-commutativity, e"*'^'^''"'"^"-' is not a unitary Anderson model and, vice 
versa, D^S can generally not be written in the form e~*(''''+^^). 

For S we choose what we consider to be a unitary analog of the discrete Laplacian. For 
d = 1 it has a five-diagonal structure which finds its roots in some physics models, see |13] . 
It turns out that S has the same five-diagonal structure as the so-called CMV-matrices, 
which have recently found much interest in the theory of orthogonal polynomials on the 
unit circle, where they arise as unitary analogs of Jacobi matrices, see [50l [52] and 
references therein. For d > 1 we define 5" as a d-fold tensor product of its one-dimensional 
version. For details see Section [2j One of the reasons for this choice of S is that one can 
view the CMV-matrix structure as the simplest non-trivial band structure which a unitary 
operator on ^^(Z) can have, see e.g. |13] . similar to the role of the discrete Laplacian 
among selfadjoint band matrices. Moreover, we choose S such that it will be invariant 
under translations by multiples of 2, yielding ergodicity of the unitary random operator 

Monodromy operators of the form (jl.2p . though not necessarily incorporating Anderson- 
type randomness, have also been proposed and studied in the physics literature [44], [6], 
[11] (see also [18]). A mathematical investigation of these models was initiated in [13] 
and continued in [lOl HH [35l |36l HZ]. In particular, spectral localization for the unitary 
Anderson model was established in [35] for the one-dimensional model and in [H] for 
arbitrary dimension in the presence of large disorder. Here spectral localization refers to 
the property that U^^ has pure point spectrum for almost every uj. 

We will provide proofs of dynamical localization (as formally defined in Section [3] be- 
low) for the unitary Anderson model in three different regimes: At arbitrary disorder and 
throughout the spectrum for the one-dimensional model as well as in the large disorder 
and band-edge regimes in arbitrary dimension (see Section [3] for a detailed description of 
these regimes). This coincides with the regimes where localization has been found to hold 
for selfadjoint Anderson models. 
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Our approach to localization proofs will be via a unitary version of the fractional mo- 
ment method, which was initiated as a tool in the theory of selfadjoint Anderson models 
by Aizenman and Molchanov in [J]. Dynamical localization will follow as a general con- 
sequence of exponential decay of spatial correlations in the fractional moments of Green's 
function (Section [5|). To complete the proof of dynamical localization in the three regimes 
described above, the latter property of Green's function will then be established in those 
regimes. 

In fact, in the large disorder regime this has already been done in ^Ij. Its proof for 
the one-dimensional model is one of the main results of the thesis |34j . from where we 
borrow the proofs presented here (Section [7]) . Some of our general results, in particular the 
proof that exponential decay of fractional moments of Green's function implies dynamical 
localization (Section [5]) and the proof that fractional moments of Green's function are 
bounded (Section [4]), are also essentially taken from [34]. 

The hardest, but possibly also most rewarding part of our work, is the proof of expo- 
nential decay of fractional moments of Green's function in the band edge regime, which is 
carried out in Sections [8] to [HI Several preparatory sections are devoted to building up 
various mathematical tools which do not seem to be known in the context of unitary oper- 
ators, such as the Feynman-Hellmann theorem from perturbation theory (Section [9]) and 
Combes-Thomas type bounds on eigenfunctions (Section llip. Along the way to localiza- 
tion we establish the spectral theoretic precursor of Lifshits tails of the integrated density 
of states for the unitary Anderson model (Section [T2|) as well as a decoupling procedure 
required in the iterative proof of exponential decay of fractional moments near the edges 
of the spectrum (Sections [TOl and [T3|) . 

In Section [6] we also include a proof of the general fact that, in the context of the unitary 
Anderson model, dynamical localization implies spectral localization, as previously known 
for selfadjoint Anderson models. In the unitary context, this follows from a version of the 
RAGE-Theorem provided in [30], whose proof we reproduce here. 

As already mentioned, when d = 1 and when restricted to /^(N) our unitary Anderson 
matrices C/^ bear close resemblance with the CMV matrices in the theory of orthogonal 
polynomials on the unit circle, see ^O]. These polynomials are determined by an infinite 
set of complex numbers on the unit disc that are called Verblunsky coefficients. Actually, 
Uuj corresponds to a choice of Verblunsky coefficients characterized by constant moduli r 
and correlated random phases, see [35] for details. Other choices of random Verblunsky 
coefficients have been studied in the literature, see e.g. [51], [52], [55] and references therein. 
We note that for i.i.d. Verblunsky coefficients in the unit disc with rotation invariant 
distribution, Simon proves dynamical localization in [51]. While not spelled out explicitly, 
our results for the one dimensional case show that dynamical localization also holds for 
the CMV matrices considered in [35] with constant moduli Verblunsky coefficients and 
correlated phases. 

Finally, there is an underlying pedagogical goal to our paper: We use the unitary 
models considered here to give a self-contained presentation of the mathematical theory of 
Anderson localization via the fractional moment approach. Making use of state-of-the-art 
techniques from localization theory, we revisit the peculiarities of the one-dimensional case 
and techniques covering various regimes in the multi-dimensional case within a unitary 
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framework. This requires developing and adapting all necessary background, which we do 
in a widely self-contained fashion. 

Acknowledgments: E. H. would like to acknowledge support through a Junior Re- 
search Fellowship at the Erwin Schrodinger Institute in Vienna, where part of this work 
was done. Also, E. H. and G. S. would like to express their gratitude for hospitality at the 
Institut Fourier of Universite de Grenoble during visits at crucial stages of this work. 



2 The Unitary Anderson Model 

As the unitary Anderson model we denote a unitary random operator of the form 

= D^S (2.1) 

in Motivated by earlier investigations in the physics literature, e.g. [llj, this model 

was studied mathematically in [13] , [35] , [IQ] , [JT] and [36] , from where we take the following 
definitions and basic results. 

A deterministic unitary operator S on /^(Z'^), sometimes referred to as the "free" unitary 
operator or "unitary Laplacian", is constructed as follows: 

Starting with d = 1, let Bi and B2 be unitary 2x2 matrices given by 



Si 



r t 
— t r 



and B2 



r —t 
t r 



(2.2) 



with the real parameters t and r linked by + = 1 to ensure unitarity. Now let Ue be 
the unitary matrix operator in /'^(Z) found as the direct sum of identical i?i-blocks with 
blocks starting at even indices. Similarly, construct Uq with identical i?2-blocks, where 
blocks start at odd indices. Define Sq = UgUo, which will serve as the operator S in (12. ip 
for dimension d = 1. The operator Sq is unitary, with band structure 
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(2.3) 



where the position of the origin in Z is fixed by {e2k-2\Se2k) = —i^-, with {k G Z) denoting 
the canonical basis vectors in /^(Z). Note also that Sq is invariant under translations by 
multiples of 2. Due to elementary unitary equivalences it will suffice to consider < t, r < 1. 
Thus So is determined by t. We shall sometimes write S'o(t) to emphasize this parameter 
dependence. The spectrum of SQ{t) is given by the arc 



a(5o(t)) = S(t) = {e*'^:T?G[-Ao,Ao]}, 



(2.4) 
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with Aq := arccos(r^ — t^). The spectrum is symmetric about the real axis and grows from 
the single point {1} for t = to the entire unit circle for t = 1. The spectrum is purely 
absolutely continuous for t > 0. 

To define the multidimensional unitary Laplacian, we follow [41] in viewing l'^(lJ^) as 
<S?'j=i^^(Z) so that for all k G Z'^, eu — e^^ ... e^^. Using Sq = So{t) from above we 
define S = S(t) by 

=0,ti5o(t). (2.5) 

The spectrum of S{t) is 

a{S{t)) = {e'^ : E [-d\o, dXo]}. (2.6) 

Throughout this paper | • | will denote the maximum norm on Z*^. Using this norm it is 
easy to see that S{t) inherits the band structure of 5*0 such that 

{ek\S{t)ei) = if|fc-/|>2. (2.7) 

For the definition of the random phase matrix in ()2.ip . introduce the probability 
space (0,9", P), where = T^'' (T = R/27rZ), 3^ is the cr-algebra generated by cylinders 
of Borel sets, and P = (S)k£Z'' ^' where /i is a non trivial probability measure on T. The 
expectation with respect to P will be denoted by E. We will assume throughout that /i is 
absolutely continuous with bounded density, 

dfi{e) = T{e) de, r e l°°(t). (2.8) 

The random variables 9k on (17, 9", P) are defined by 



7fc 



In other words, the {0%}ke'L''- T-valued i.i.d. random variables with common distribution 

/X. 

The diagonal operator in £^(Z'^) is given by 

D^ek = e-'^^efc. (2.10) 
With this choice for and S = S{t) we define the unitary Anderson model C/o; via 

This definition and the periodicity of S ensures that the operator is ergodic with 
respect to the 2-shift in fi. Uuj also inherits the band structure of the original operator 
S. The general theory of ergodic operators, as for example presented in Chapter V of [17j 
for the self-adjoint case, carries over to the unitary setting. In particular, it follows that 
the spectrum of U^j is almost surely deterministic, i.e. there is a subset S of the unit circle 
such that a{Ui^) = S for almost every w. The same holds for the absolutely continuous, 
singular continuous and pure point parts of the spectrum: There are 'Sac, '^sc and T,pp such 
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that almost surely (Tac{Uu)) = Sqc, (^sc{Uuj) = ^sc and crpp{U^) = Spp. Moreover, we can 
characterize S in terms of the support of /j, and of the spectrum of 5; 

S = exp {—i supp fi) a{S) = {e*° : a £ [—d\o, d\o] + supp/i}. (2-11) 

Here supp /i denotes the support of the probability measure fi, defined as 

supp^ := {a I fj,{{a - e, a + e)) > for all e > 0}. (2-12) 

The identity (|2.1ip is shown in [40j for the one-dimensional model, but the argument carries 
over to arbitrary dimension. 

As i ^ 0, 5*0 (t) tends to the identity operator, whereas as f — > 1, So{t) tends to a direct 
sum of shift operators. Accordingly, if t is zero then the operators Ui^ are diagonal and 
thus trivially have pure point spectrum. On the other hand, if t = 1 then it is not hard to 
see that all [/^ are purely absolutely continuous (in fact, they are unitarily equivalent to a 
direct sum of shift operators). Thus, excluding the trivial special cases, we shall from now 
on assume that < t < 1. For the unitary Anderson model the parameter t takes the role 
of a disorder parameter with small t corresponding to large disorder, as this means that 
C/^ is dominated by its diagonal part. 



3 The Results 

As discussed in the introduction, our main goal is to study regimes in which a quantum 
mechanical system governed by the unitary propagator Ui^ is dynamically localized. This 
can be expressed in terms of the transition amplitudes {ek\U^ei), whose squares measure 
the probability that a system initially in state e; evolves into state after time n. By 
dynamical localization we will refer to the property that the expectation of these amplitudes 
stays exponentially small in the distance of k and /, uniformly for all times, i.e. the existence 
of constants C < oo and a > such that 

E(supKefc|f/>z)|)<Ce-"l'^-'l. (3.1) 

In fact, what we will prove in several corresponding regimes is the stronger result that 

E(supKefc|/([/^)Q)|)<Ce""l'=-^ (3.2) 
/ 

where the supremum is taken over all functions / G C{S) with ||/||oo := sup^gs \ fiz)\ < 1. 
Here S = {z G C : |z| = 1} is the unit circle. 

Also, dynamical localization may only hold in an arc {e*^ : £ [^)^]} of the spectrum 
of Uuj- In this case, the spectral projection -Pj-'^^j of onto this arc will be applied to the 
state ei in (j3.2p . restricting the initial state to the localized part of the spectrum. 

Our detailed results, stated in the following two subsections, fall into two categories: 
We start with results which show that dynamical localization can be established by a 
unitary version of the fractional moments method, using as a criterion the exponential 
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decay of fractional moments of spatial correlations of the Green function. We also show 
that dynamical localization generally implies spectral localization, i.e. that U^j almost surely 
has pure point spectrum in the corresponding part of the spectrum. Finally, we state that 
dynamical localization implies almost sure finiteness of all quantum moments of the position 
operator. All these results hold for arbitrary dimension d, for general values of the disorder 
parameter t E (0,1), and without restriction on the spectral parameter of the unitary 
operators U^. 

Our second set of results concerns the proof of dynamical localization in three concrete 
regimes: for the one-dimensional unitary Anderson model, as well as for large disorder 
and at band edges in arbitrary dimensions. In each case this will be done by verifying 
exponential decay of the fractional moments of the Green function. 



3.1 Fractional Moment Criteria for Localization 

For z e C with \z\ / 1 let 

G{z)=G^{z) = {U^-z)-\ (3.3) 

and 

G{k, I] z) = {ek\G{z)ei), k,leZ'^ (3.4) 

be the Green function of (to use a term from the theory of selfadjoint hamiltonians in 
the unitary setting). 

The Green function becomes singular as z approaches the spectrum of Ui^. The first 
insight which makes the fractional moment method a useful tool in localization theory is 
that these singularities are fractionally integrable with respect to the random parameters. 
This means that for s G (0,1) the fractional moments K{\G{k,l; z)\'^) have bounds which 
are uniform for z arbitrarily close to the spectrum. This is the content of our first result. In 
fact, we will need a somewhat stronger result later, namely that it suffices to average over 
the random variables Of. and 9i to get uniform bounds on \G{k, I; z)\^. The role of bounds 
of this type in localization proofs via the fractional moment method roughly corresponds 
to the use of Wegner estimates in the approach to localization by the method of multi scale 
analysis. 

Theorem 3.1. Assume that the random variables {9k}k&'^ '^'^^ i.i.d. with distribution ^ 
satisfying ^2.3\) . Then for every s G (0, 1) there exists C(s) < oo such that 



\G{k,l-z)Ydii{ek)dii{9i)<C{s) (3.5) 

for all z G C, \z\ ^ 1, all k,l G Z"', and arbitrary values of 9j, j {k,l}. Consequently, 

Ei\Gik,l;z)n<Cis), (3.6) 
for all z £ C, \z\ ^ 1. 
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The second statement simply derives from the bound (j3.5p . which uniform in the random 
variables 9j, j {k,l}, and the independence of the 6j. The proof of (j3.5|) is given in 
Section [H 

In the next subsection we will identify several situations where the fractional moments 
K{\G{k, I; 2)1*) are not just uniformly bounded, but decay exponentially in the distance of k 
and I. The following general result shows that this can be used as a criterion for dynamical 
localization of Ui_j. 

Theorem 3.2. Assume that the random variables {0k}keZ'' satisfy 1^2. 8\) and that for some 
s G (0, 1), C < 00, a > 0, e > and an interval [a, b] G T, 

E(|G(A:,/;z)n < Ce-"l*--'l (3.7) 

for all k,l 7j'^ and all z ^ C such that 1 — e < |2;| < 1 and argz G [a, b]. 
Then there exists C such that 

E[ sup |(efc|/(C/.)P[-,]e,)|]<Ce--l^-^l/^ (3.8) 
/ec(s) 

ll/l|oo<l 

for all k,l£Z'^. 

Here, as usual argz G T refers to the polar representation z = |z| exp (i arg z) of a 
complex number. We prove Theorem 13.21 in Section [5l 

It has been shown in ^41j that fractional moment bounds of the form (j3.7p for a unitary 
Anderson model imply spectral localization via a spectral averaging technique, following an 
approach to localization which is due to Simon and Wolff [5l] for the selfadjoint Anderson 
model. For completeness, we present a direct proof of the fact that dynamical localization 
expressed by ()3.8p implies spectral localization in the unitary setup. This follows from 
a simple adaptation of arguments borrowed from Enss-Veselic |30j , see also [15], on the 
geometric characterization of bound states, i.e. a RAGE-type theorem for unitary operators. 
We prove 

Proposition 3.1. Assume that for an interval [a,b] there exist constants C < 00 and 
a > such that 

E[ sup Kefc|/(C/jP[-,]e,)|]<Ce-"l^-^l (3.9) 
/ec(s) 

ll/l|oo<l 

for allkJeZ"^. Then 

(a, 6) n Scont = 0, (3.10) 
where Ticont = U T^ac- In other words, almost surely Py^^^^U^j has pure point spectrum. 
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Another direct consequence of the dynamical locahzation estimate ()3.8p . is that it pre- 
vents the spreading of wave packets from P\^^ under the discrete dynamics generated 
by Uu)- This dynamical localization property is measured in terms of the boundedness in 
time of all quantum moments of the position operator on the lattice. More precisely, for 
p > let |X|e be the maximal multiplication operator such that 

\X\lej = \j\lej, forjGZ^ (3.11) 

where \j\e denotes the Euclidean norm on 7J^. We have 

Proposition 3.2. Assume that there exist C < oo and a > such that 

E[ sup |(efc|/(C/jP[-,]ei)|]<Ce-l^-^l (3.12) 
/ec(§) 

ll/l|oo<l 

for all A;, / G Z'^. Then, for any p> and for any ijj in P{'Z'^) of compact support, 

sup|||X|Pf7^P[';^ ^Vll < oo almost surely. (3.13) 

Similar results hold under weaker support conditions on ip. Our choice is made to keep 
things simple. The proofs of Propositions 13.11 and 13.21 are given in Section [6l 

3.2 Localization Regimes 

The other main goal of our work is to identify three different regimes, where the frac- 
tional moment condition ()3.7p can be verified, and thus dynamical localization follows by 
Theorem [321 



3.2.1 Large Disorder 

As explained at the end of Section [2l the parameter t E (0, 1) used in the definition of 
S{t) and thus, implicitly, also in the definition of Ui^, can be thought of as a measure of 
the degree of disorder in the unitary Anderson model. Thus one can expect a tendency 
toward localization for small values of t. This was confirmed in [51] where the following 
was proven: 

Theorem 3.3. Suppose that the i.i.d. random variables {0k}keZ'' have distribution jj, sat- 
isfying \2. (S'p and let s E (0, 1). Then there exists to > and C < oo such that if t < to, 
there exists a > so that 

E{\{ej\UUUu. - z)-'ek)n < Ce-"^'-"^ (3.14) 
for all j, k and all z with \z\ / 1. 
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In fact, [41j considers a more general model in which different parameters ti are chosen 
in each factor of (12. Sp and shows (|3.14p under the condition that U is sufficiently small. 

Using (j4.ip below this implies that (j3.7p holds for all \z\ ^ 1. Thus in the large disorder 
regime t < to dynamical localization holds on the entire spectrum of U^j by Theorem 13.21 
(the spectral projection Pr'^y in (|3.8p can be dropped). 



3.2.2 The One-Dimensional Model 

For the one-dimensional self-adjoint Anderson model, localization holds throughout the 
spectrum, independent of the amount of disorder. The same is true for the unitary Ander- 
son model, as implied by the following result: 

Theorem 3.4. Let d = 1 and suppose that the i.i.d. random variables {9k}kGZ have distri- 
bution fj, satisfying \2.^) . Then for every t < 1 there exist s > 0, C < oo and a > such 
that 

E(|G(A:,/;z)n < Ce-"l^-^l (3.15) 
for all z e C such that < ||2;| - 1| < 1/2 and all k,l G Z'^. 

By Theorem 13.21 this implies dynamical localization for the one-dimensional unitary 
Anderson model throughout the spectrum. 

Many of the special tools which have been heavily exploited in studies of the one- 
dimensional self-adjoint Anderson model, are also available for the one-dimensional Unitary 
model. First of all, there is a transfer-matrix formalism which allows the definition of 
Lyapunov exponents. In particular, it has been shown in [36] that under assumption (12. 8p 
(in fact, for much more general distributions) the Lyapunov exponent is positive on the 
entire spectrum and continuous in the spectral parameter. This is the central ingredient 
into Lemma 17.11 stated in Section [7] below. For a proof of Lemma 17. H which is a close 
analog of a result proven for the selfadjoint one-dimensional Anderson model in [16], we 
will refer to [3lj . In Section [7] we will then explain in detail how this leads to (j3.15p . 



3.2.3 Band Edge Localization 

For notational simplicity and without restriction we assume for the following result that 
supp/i C [—a, a] with a G (0, vr) and —a, a G supp//. Furthermore, we assume that 

a + dXo<TT, (3.16) 
which by (j2.1ip guarantees the existence of a gap in the almost sure spectrum S of Ui_j, 

{e'^ : 1? G (dAo + a, 2^ - dAo - a)} n S = 0, (3.17) 
and that e*('^'^o+") and Q^i'^'^-'i>'o~a) band edges of S. Our main result is 
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Theorem 3.5. Assume \2. 8\) and let < s < 1/3. There exist 5 > 0, a > and C < oo 
such that 

E(|G(A:,/;z)n < Ce-"l^-^l (3.18) 

for all k,l G Z'^ and all z £ C with \z\ ^ 1 and argz G [JAq + a — 5, dX^ + a] U [27r — JAq — 
a, 27r — dAo — a + S\. 

Note that by Theorem 3.2 this imphes dynamical localization for near the edges 
dAo + a and 27r — dX^ — a of its almost sure spectrum. 

The strategy of the proof of Theorem 13.51 is the following. We control the expectation 
value of fractional moments of the infinite volume Green function in terms of the expecta- 
tion value of fractional moments of the Green function of a finite volume restriction of the 
operator. This requires addressing several distinct issues. 

The first one is the definition of an appropriate finite volume restriction. We restrict 
the problem to a finite but large box, by introducing appropriate boundary conditions 
in Section 8. Our choice of boundary conditions is governed by the fact that we need 
monotoncity properties which are similar to those of Neumann conditions in the selfadjoint 
case. The link between the infinite and finite volume resolvents is provided by a geometric 
resolvent estimate and a decoupling argument, similar to the self-adjoint case. Provided 
one has good estimates on the expectation value of the fractional moments of the finite 
volume Green function, this allows to lift such estimates to the fractional moments of the 
infinite volume resolvent by means of an iteration, for large but fixed size of the box, in a 
neighborhood of the band edges. This second step is addressed in Section [I3l 

To get the sought for estimates on the resolvent of the finite volume restriction, in the 
band edge regime, we need to control the probability that the finite volume restriction of 

has eigenvalues close to the band edge. Quantitatively, this amounts to showing that 
the probability of a small distance, algebraic in the inverse size of the box, between the 
eigenvalues of this finite volume restriction and the band edges is exponentially small, as 
the size of the box increases, see Proposition [T27TJ This is an expression of the fact that the 
spectrum close to the band edges is very fluctuating which gives rise to Lifshits tails in the 
density of states. To prove this, we follow the self-adjoint route, see [56j. We first study the 
effect of Neumann boundary conditions on the spectrum of S in Section [HI and we make 
use of a unitary version of the Feynman-Hellmann formula. Proposition [9Tl to control this 
effect on the spectrum of the random operator Ui^. Then, the Lifshits tail estimate together 
with a unitary version of the Combes-Thomas estimate, Proposition 111.11 allow to show 
that the expectation of the moments of the finite volume Green function is exponentially 
small in a power of the size of the box. Proposition 114.11 

4 Boundedness of Fractional Moments 

In this section we prove Theorem 13.11 
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As the bound (j3.5p is trivial for \z\ < 1/2, it suffices to consider \z\ > 1/2 and \z\ ^ 1. 
Below we prefer to work with the modified resolvent {Uuj + z){Uuj — z)""^ . Since 

{U^ - z)-^ = ^[([/^ + z){U^ - z)-^ - I], (4.1) 
it easy to see that the existence of < C < cxd for which 

\{ek\{u^ + z){u^-z)-^el)\'d^Ji{ek)d^i{el)<c, (4.2) 



for all z with \z\ 7^ 1, all k, I G Z"^, and uniformly in 9j, j {k, I}, gives the required bound. 

Key to the proof of (|4.2p is knowledge of the exact algebraic dependence of the Green 
function on the two parameters 9^ and 61. Similar formulas for rank one and rank two 
perturbations of the resolvents of unitary operators have been derived in |49j, Section 4.5, 
from where we took guidance. 

We mostly focus on the proof of (|4.2p for the case k ^ I. At the end of the proof we 
comment on the simpler case k = I, where ()4.2p only requires averaging over one parameter. 

For k 7^ / weborrow an idea from |3j and introduce the change of variables a = ^{9k+9i), 
/3 = ^(^fc ~ (^i)- This will have the effect of essentially reducing (j4.2p to averaging over the 
single parameter a (although this still corresponds to a rank two perturbation). 

Let A = {k, I] C Z'^ and define 

= { -f3, j = l (4.4) 

(4.5) 

Next, we define the diagonal operators D^, Dp and D by 

DaCj = e'^'^^ej, DpCj = e'^^^ej, Dej = e^^'^^Cj. (4.6) 

Using these definitions we can write 

= D^V^, (4.7) 

with the unitary operator = D^DS. In what follows we explore the relation between 
the modified resolvents of U^j and V^^. 

Let Pa be the orthogonal projection onto the span of {V^^ej : j G A}. Using that 
{V~^ej : j G Z'^} is an orthonormal basis of /^(Z*^), simple calculations show that (C^j — 
K))(/ — Pa) = and Vj^Ui^ = e~*"I on range Pa- In particular, — Vuj = {U^j — K;)-Pa 
is a finite rank operator. Therefore, 

= VM - Pa) + e-'"K.PA. (4.8) 
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For z£C\{0} with \z\ + 1, let = P AiJJ ^ z){\J ^ - z)-^ P a while F^, = + 
z)^^Pa^ both viewed as operators on the range of Pa (i-e. 2 x 2-matrices). We see that 

% + F; = Pa(2/ - 2|zp)(K. - ^)~H(K. - zY^XPa- (4.9) 

This shows that F^ + F* is invertible and + F* < for > 1. Therefore, —iF^ is a 
dissipative operator, i.e. an operator A such that (A — A*)/2i > 0. Similarly, —iF~^ is 
also a dissipative operator. In the case |2;| < 1, we have that iF^, iF~^ are dissipative. 

Next we explore the relation between F^ and F^. Following Section 4.5 of [39] we use 
the fact that {x + z){x — z)~^ = 1 + 2z[x — z)^^ along with (14. Sh to obtain 

F,-F, = -2zPa{V^ - zy^VM^'"'" - 1)Pa{U^ - zy^PA. (4.10) 
Using the definitions of F^ and F^, this can be rewritten in the form 

F,-F,= '^{1 + F.)(e-'- - 1)(1 - F,). (4.11) 

For a ^ {0,7r}, let m(a) = i |^j]g-iQ G M. A straightforward calculation shows that 

F^ = -i{-iF^ + m(a))"^ - ii-iF'^ - m"^(a))"^ (4.12) 

Note that F^ depends on (3 and the 9j, but not on a. 

From the definitions of F^ and Pa and the fact that V~^Uij = e~*"I on span{T/7^ej : 
j G {A;, Z}}, a simple calculation shows that 

{ek\{U^ + z){U^-z)-^ei) = {V-^ek\F,V-^ei). (4.13) 

Therefore, 

I (efc|(C/^ + z){U^ - zy\i)Y dfiiOk) dii{9i) 

'2n f2n 



< \\r\\lo r \ {V-^ek\F,V^^ei)\' dOkdOi 

Jo Jo 

/'27r /'27r 

< \\r\\lo / / llF^rd^fcrf^z 



r-2n 

|2 



'0 JO 

< 2||r||^ lli^.rdaj (4.14) 

where we have changed to the variables a and (3 and slightly enlarged the integration 
domain into the rectangle < a < 2tt, —tt<P<tt. 
We split the a-integral according to ()4.12p . 

f27T p2iv /•27r 

/ \\Fz\Yda< \\{-iF^ + m{a)y^\Y da+ ||(-iF-^-m-^(a))-^f da. (4.15) 
Jo Jo Jo 
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Recalling that m(a) has singularities at a G {0,27r}, we make the change of variables 
X = m(a), 



l'2-iT 

Jo 



-iFz + m{a)) \[''da 



lim / -J—\\(-iF, + x)-H\'dx 



R~*oo 7„^ + 1 
n+1 



2E / ^iiHi^.+-)-^irdx 



We can now complete the proof of (j4.2p by treating the cases |2;| > 1 and |2;| < 1 separately. 
If |z| > 1 then —iFz is dissipative and Lemma l4.ll shows boundedness of the integral on 
the right of (j4.16p . uniform in n, \z\ > 1, (3 and dj. Thus, after summation, j"^^ \\{—iFz + 
m(a))~^||* < C{s). The second term on the right of (j4.15p can be bounded in a similar 
way, using that —iF~^ is dissipative as well. Inserting these bounds into (j4.14p makes the 
/3- integration trivial and completes the proof of ()4.2p for the case \z\ > 1. 

li\z\ < 1, then —iFz and —iF^^ and anti-dissipative. As Lemma [4.11 obviouslv also holds 
for anti-dissipative matrices j4, the proof of (14. 2p goes through with the same argument. 



The proof of (j4.2p for the case k = I is similar but simpler. We don't need a change of 
variables, but directly work with one of the parameters 9i instead of a, leading to rank one 
perturbations. The objects corresponding to F^ and Fz become scalars and we only have 
to use the trivial scalar version of Lemma |4. II to conclude. 

We finally provide an elementary proof of the following Lemma which was used above. 
A much more general result of this form is given as Lemma 3.1 in j3]. 

Lemma 4.1. For every s £ (0, 1) there exists C(s) < oo such that 

\\{A + xI)-^\\' dx <C{s) (4.17) 

E 

for every dissipative 2 x 2-matrix A and every unit interval E. 

Proof: First observe that a general dissipative 2 x 2-matrix A is unitarily equivalent to 
an upper triangular dissipative matrix (choose as unitary transformation any matrix whose 
first column is given by a normalized eigenvector of A) . Thus we may assume that 

which implies 

/ 1 ai2 \ 

{A + xI)-^=[ ""+^ (aii+x)(a22+x) \ ^419^ 

V 122+3; / 
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The bound ()4.17p follows if we can establish a corresponding fractional integral 
for the absolute value of each entry of (14.19P separately. For the diagonal entries 
obvious. 

We bound the upper right entry of (j4.19p by 



bound 
this is 



ai2 



(an + x)(a22 + x) 



The positive matrix 

lmA = ^{A- A*) 



< 



ai2 



|Im((aii +x)(a22 + x))\ 



. Imaii+Imc 
|ai2| 



Imaii Jjai2 



+ 



Im(aiia2 



2^012 Ima22 



has positive determinant, i.e. det Imj4 = Imaiilma22 — |ai2|^/4. We thus get 

2 



Imaii + Ima22 



0.12 



^ 21maiilma22 ^ 1 



(4.20) 



(4.21) 



(4.22) 



The latter allows to conclude the required integral bound for (j4.20p . 



5 Dynamical Localization via Green's Function 

Here we will prove Theorem l3.21 i.e. that exponential decay of fractional moments of Green's 
function implies dynamical localization. 

Our proof uses an idea which in the context of selfadjoint Anderson models is due 
to Graf [33], namely that second moments of an Anderson model's Green function can 
be bounded in terms of its fractional moments (including, however, a scalar factor which 
becomes singular as the spectral parameter approaches the spectrum). While the details 
of the proof are more involved than in the selfadjoint case, we find a bound of this form 
for unitary Anderson models in Section 15. 1[ 

Another tool we use is the integral formula (j5.28p . which expresses operator functions 
/([/) in terms of the resolvent of U . In Section 15.21 we provide a proof of this formula, 
which combines the spectral theorem for unitary operators with the representation of Borel 
measures on T by Poisson integrals. 

Equipped with these tools we complete the proof of dynamical localization in Sec- 
tion (531 

5.1 A Second Moment Estimate 

We start by a bound of second moments of Green's function in terms of its fractional 
moments, which holds pointwise in the spectral parameter. 
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Proposition 5.1. Assume that the {^^ j^ez'^ satisfy i2.8]) . Then for every s £ (0, 1) there 
exists C{s) < oo such that 

Ei{l-\z\^)\Gik,l;z)\')<Cis) n\Gim,l;z)n (5.1) 

|m-fc|<4 

for all \z\ < 1 and k, I € Z"^. 



Proof: Throughout the proof we will assume z ^ 0. The bound (j5.ip carries over to 
z = by continuity. 

For 6 eT, let % = e-'^^k+^) - e'^^k . Then define = + r]kPk, where Pk is the 
orthogonal projection into the span of e^. Let = D^S. Using the resolvent identity we 
have 

{U^ - z)-^ - {U^ - = -{Ut - z)-%PkS{U^ - z)-\ (5.2) 

for all z G C such that < |z| < 1. Letting F{z) = S{U^ - z)-^ and Fs{z) = S{U^ - z)-^ 
the last equation takes the form 

Fs{z)-F{z) = -Fs{z)r]kPkF{z). (5.3) 

Denoting F{i,j,z) = {ei\F{z)ej) it is easy to see that 



Therefore, for all / G 



7d 



On the other hand, we also have that 

\Fs{k,l,z)\^ < Y\Fs{k,y,z)\'' 

yez 

= {ek\SiU^-z)-\iU^-z)-'rS*e,). (5.6) 
Since is a unitary operator, the following identity holds 

[{U^-z)-'r = =^iU^-lj-'u!,. (5.7) 
Thus, it follows that 

\Fsik,l,z)\^ < =l{ek\S{Ut-zr\U^-h-'Dtek) 



z z 
_e-i(0^+5) 
— {ek\S{U,:,- z) ^U:,- 

z 



{ek\S{U^-zr\U^--)-'ek). (5.8) 
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Again using the resolvent identity, we see that 

(U^ - zr\ut - \r' = T^^iiU^ - z)-' - {Ut - \)-'}. (5.9) 



Hence, 



-j(e^+(5) 1 

\F5{k,l,z)\'' < {{ek\S{U^ - z)-'ek) - {ek\SiU^ - ^y'ek)}. (5.10) 

1 — zr z 



Prom the definition of and the fact that (U^ - z)-^ = -^[I - U^{U^ - z)-^], it 

follows that 

{ek\Siu!,-^)-'e,) = -ze'(9t+^Hek\mt-z)-']*ek) 

= -ze-(^^+^){l-(^^ ([/i-z)-^ efc|efc)} 

= e^(^^+'5){l-e^(''^+^)F5(A;,A;,z)}. (5.11) 

Therefore, we now obtain that 



mk,l,z)\^ < -^-^—{2?fi[e-'('"^+'^Fs{k,k,z)]-l} 



= -r^i\Ps{k,k,z)\^ - \e''^'^^+'^ -Fs{k,k,z)\^}, (5.12) 
1 — |z|^ 

since \x — = |xp + |yp — 23f?[2;y]. Using (j5.5p . to rewrite Fs{k, k, z) in terms of elements 
of F, along with the definition of tj^ we get 

lj,fl.l ^,2^ 1 , \Fik,k,zr-\e^'^^ -F{k,k,z)\\ 

l'^^^^'''"^' -T^^ |l + ,,F(fc,fc,z)P ^'-''^ 

This inequality gives, in particular, that F{k, k, z) ^ 0. Therefore 

Finally note that the last inequality allows us to conclude that |1 — e^^''F{k, k, z)~^\ < 1. 
One can also use the fact that 

\Fs{k,k,z)\ < \\{U^ - z)-i| < (5.15) 

for all (5 G T, to get a different upper bound on (1 — \z\'^)\Fs{k,l, z)\'^ . Since (|5.5p can be 
rewritten as 

' ' ^ rjk + F(k,k,z)-^ F{k,k,z)' ^ ' 

it follows that 1 — |% + F{k,k, z)^^\ < \z\. Then by choosing 6 such that e"**^ = 
T- — »o, c^;, , x_i , we see that 

\l-e''^^F{k,k,zy^\ <\z\. (5.17) 



18 



Using this along with (j5.16p we obtain the following upper bound 



{l-\z\')\Fs{Kl,z)\' < 



2 e''^''F{k,k,z)-^\^ \F{k,l,z 



|% + F(A;,A;,z)-i|2 \F{k,k,z)\^' 



(5.18) 



Combining the two estimates ()5.14p and (|5.18p and using that for < s < 1 we have 
min(l, < it follows that 



(l-|z|2)|F5(A;,/,z)|2 < 



1 - |1 - e'^''F{k,k,z) 



-1|2 



\F{k,l,z)Y 



-iS 



(l-e*«fcF(fe,fc,z)-i)|2 \F{k,k,z) 



Letting y = 1 — e^^''F{k, k, z) ^, this can be rewritten as 

(1 - \z\'mik,i,zr < ^']^fsl\y^^ \Fik,l,z)\^ 

Since the expectations of F and Fg are related by 
E[\F{k,l,z)\^]=E[[ dfi{9k + 6)\Fs{k,l,z)W 



(5.19) 



(5.20) 



(5.21) 



it follows that 

E[il-\z\^)\Fik,l,z)\^ 

< ||r||ooE 



\F{k,l,z)\' sup 

{y£C:\y\<l}JO 



2-^^(i-M')|i-yr 



-iS 



< 2l|r|UE 



\F{k,l,z)\' sup 

{yeC:\y\<l}J0 



2tt 



d5 



(l-bp) - 



(5.22) 



Next we evaluate the integral 



-iS 



2/P 



27r -iS I 

d6^ 







27r 



K / (i(5[ 







g iS _ y 

. 26-*"^ 

g 10 _ y 



!]■ 



(5.23) 



The latter integral can be easily evaluated using Cauchy integral formula, by simply sub- 
stituting z = e~^^ and gives 

27T. (5.24) 



3— i(5 



(5.25) 



Therefore, 

E[(l - \z\^)\F{k,l,z)\^] < r+'7r\\T\\^E[\F{k,l,z)n 
Since S" is a unitary operator with band structure, we see that for s S (0, 1) 
E[\F{k,l,z)n = E[\{S*ek\{U^-zr'ei)n 

< Ci{s) mem\{U^-z)-'ei)n. (5.26) 

\m~k\<2 
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Finally, in order to get the required bound of the second moment of elements of (Uuj — 
z)~^ we use again that S" is a unitary operator with band structure. Therefore, for k,l £ 1^'^ 

E[{l-\z\^)\{ek\(,U^-zy'ei)\^] = E[{1 - \zf)\{Sek\S{U^ - zr'ei)f] 

< C2is) Yl m-\z\')\Fim,l,z)\\ (5.27) 

\m-k\<2 

Combining (fOTj) . (fOSll and (lOGll gives (^J^. ■ 



5.2 An Integral Representation for f{U) 

We will reduce bounds for f{U) to bounds on resolvents by the formula 



f{U) = w- lim / {U -re''')-\V -re-'')-'f{e'')de (5.28) 

r^l- 27r Jo 

for / E C(S) and U a unitary operator. This is a simple consequence of the representation 
of Borel measures on T by Poisson integrals: 

Let (/3 G J{ be normalized and consider the non negative spectral measure d/i^ on T 
such that 

{ip\U^) = I e^"d(99|S(a)99) = / e^"d/i^(a), (5.29) 
where E(-) denotes the spectral family of U. We can thus rewrite for < r < 1 



I 

JT 



1 -r^ 



1 + — 2r cos{9 — a) 



1 - 

J |e'" — re" 

dfj.^{a) = u{r,6), (5.30) 



which coincides with the Poisson integral of the measure df^^p . As a function oi z = x + iy = 
re*^, with the standard abuses of notations, the RHS of (j5.30p is non negative and harmonic 
in the open unit disc ([17] Sect. 11.17), but it is not bounded as r ^ 1~ at the atoms of 
d/j,^. Now, Theorem 3.9.8 in [l6] on the representation of non negative Borel measures on 
T says that for any / G C(S) we have 

r2-iv in !• 

lim / u{r,e)f{e'')—= f{endfi^{a)^{^\f{U)^). (5.31) 

Considering the matrix element • ip) for arbitrary (p and ip in !K, we get the equivalent 
result by polarization, which proves (I5.28p . 

If P[afi] denotes the spectral projector of U on the arc [a, 6] C T and 99 G !K is normalized, 
then 

d{^\Pya,b]E{a)ip) = a!/u[,"'^](a) = di^^{a)\ya,b]. (5.32) 
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By the same argument with -P[a,fe]V in place of y?, i.e. with di^^'^"^ in place of d/i^, we 
have for any f G CiS) 

P[a,b]fiU) = w- YmiJ^ j\u-r^')-\U-^-re''')-^f{e'')de 

= u;- lim fr{U). (5.33) 

5.3 Conclusion of the Proof of Theorem 13.21 

Let f £ C{§) such that ||/||oo < 1 and consider 

f^(U) = [\u - re^^)-i(C/-i - re-'^)-^f{e'^)de. (5.34) 

We compute 

\{ek\fr{U)e,)\ 

^ •{eu\{U -re'')-\U-^ -re-'')-^e,)\\f{e'')\de 
^ ^ r^Ke,|([/-re^^)-ie,)|Ke,-|(C/-re^^)-ie,)|d0, (5.35) 

which is independent of /. Hence, using Fatou's Lemma and Fubini's Theorem and taking 
the supremum over all / G C(S) such that ||/||oo < 1, 

E[swp\{ek\Pia,b]fiUu)ej)\] 
= E[sup lim \{ek\friU)ej)\] 

< E[liminf sup \{ek\fr(.U)ej)\] 

'■^i" ll/IU<i 

< liminfi^ /" y2E[\{ek\{U-re'')-'ei)\\{e,\{U-re''r'ei)\]\d9. (5.36) 



By the Cauchy-Schwarz inequality and the second moment bound (j5.ip . we have 
E[|(l - r2)(efc|(C/ - re^'r\i)\\{e,\{JJ - re^^)"^e,)|] 

< (E[|(l - r2)(efc|(C/ - re*^)-ie,)p])V2(E[|(l _ r2)(e,|(C/ - re*^)-iei) p]) V2 

< cA IE(G(m,Z;re*^)r) J] E(G(m, /; re*^)|^ 

\m:|m— fc|<4 / \ m:|m— j |<4 

< Cae-^l'^-'l/^e-^IJ-'l/^, (5.37) 

where ultimately the assumption (j3.7p on exponential decay of fractional moments of 
Green's function was used. Since there exists a finite c such that 
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we finally get 

E[sup\{ek\Pia,b]fiUu)ej)\] < Ce'^\^-^\l\ (5.39) 
with C = c C2 . This completes the proof of Theorem 13. 2[ 

6 Dynamical Localization implies Spectral Local- 
ization 

We start by recalling the geometrical concepts and results introduced in |30j in a general 
framework. 

Let [/ be a unitary operator in a separable Hilbert space and let J{pp(C/), respectively 
!Kc(C/) denote the pure point, respectively continuous, spectral subspace of \J . We will 
denote the spectral resolution of U by E\j . 

A practical characterization of the vectors of 3fpp(f7) and !Kc(C/) in terms of their be- 
haviour under the discrete dynamics with respect to families of finite dimensional subspaces 
of J{ is provided in pO] and reads as follows: 

Let P = {Pr}r>o be a family of projections on 'K such that 

P,, = = p*^ Rank Pr < 00, and s - lim P,. = I. (6.1) 

r— >oo 

We define the set of bounded trajectories with respect to the family P by 

M^(P)={^G?{| lim sup||(/-P,.);7"^|| = 0}. (6.2) 

Similarly, the set of propagating trajectories with respect to P is defined by 

1 ^ 

KiP) ={i^&^\ lim — — WPrU'^'^W = 0, Vr > 0}, (6.3) 

n=-N 

Both M^(P) and M^(P) are easily seen to be closed subspaces of "K. 
Then the following holds. 

Lemma 6.1. For any unitary operator U on a separable Hilbert space "K, and any family 
P as in \6.1\) . 



M^(P) = 'Kpp{U) and M^(P) = 'Kc{U) (6.4) 

While the definition above suffices for our purpose, more general families of operators 
than P can be considered, see [SO]. For completeness, we provide a detailed proof of this 
Lemma, which is a slight adaptation of the argument in [30] . 

Proof of Lemma 16.11 in three steps: 
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First, 3ipp{U) C M^(P): 
If Uip = e'"(^, then 

- Pr)U''ip\\ = - PrM\ ^ Q, uniformly inn as r ^ oo. (6.5) 



We conclude using that !Hpp{U) = {(p \ U(p = e^'^ip} and the fact that M^(P) is closed. 
Second, M.\j{P) ± M^(P): 
Let if e M^{P) and ip G M^(P). Then 

iV 1 ^ 

(^1^) = 2ivTi ^ ^^'^^ = 2ivTi ^ 

n=-N n=-N 



1 ^ 

— ^ 5; ([/>!(/ - P.)C/» + (?7>|P.i7». (6.6) 



2N + ^ 

n=-N 

By our choice of ip, for any e > there exists r(e), uniform in n, such that |(t/"(^|(/ — 
Pr{e))U""ip)\ < e- And our choice of (p and the selfadjointness of Pj.(f) imply 

TV 1 ^ 

5^ (C/>|P,(,)C/'^^) <2]vTT ^ r.(.)t^>llll^ll-0 asAT^oo, (6.7) 

n=-N n=—N 

which shows that {p\tp) = 0. 
Third, JiciU) C M^(P): 

Let V' € ^c{U) = Pc{U)J{, where Pc{U) is the spectral projector of U on the continuous 
spectral subspace of U. Since the orthogonal projector P^ is finite dimensional for any r, 
there exists vectors (pi G Ran Pj. such that we can write for some m < oo 



Pr = '^\Vi){^i\, where {ipi\^Pj) = 6ij. (6.8) 

i=l 

Hence, WPrWipf = Y^^L^ KipilWip)]^ . Now, by the Spectral Theorem we have 

/^,,^(n) := (<^.|^» = [ e""d;.^^,^,(x), (6.9) 



where d^^p^^tpix) = d{Pc{U)ipi\Eu(x)Pc{U)ip) is a continuous complex valued measure. 
Thus, by Wiener's Theorem, 



1 ^ 

iv^^oo2ivTT ^ l^»,V'WI' = El/^v^».^(W)l' = o. (6.10) 



n=-iV xeT 

Consequently, making use of Cauchy-Schwarz, 

. 1/2 
\PrU-i,f I 



AT r N 

+ 1 II -^11 - ] 2Ar + l 

n=—N K n=—I\ 

{m ^ N ^ 

i=l n=-N ) 
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which tends to zero as ^ 0. 

Finally the decomposition "K = ^pp(U) © 3ic(f/) leads to the conclusion. H 

Proof of Proposition 13.11 

Coming back to the random situation at hand, we construct a suitable family of pro- 
jectors P = {Pr}r>o by means of 

Pr= Yl '^>0, (6.12) 

|fc|<r 

and consider the vectors of the form Pj*^ ^-^ej, j G Z'^. We have for all n G Z 

IIU-^.)C^:^M]e,|| ={E l(efc|P[^:,,]Ce,)p}V2. (6.13) 

\k\>r 

Therefore, since \{ek\Pl'^^h]U>j)\ < 1, 

sup||(/-P.)C/:Pf-,je,|| = {sup5]Ke,|P[-,][/:e,)|2}i/2 

n n , 

\k\>r 

< {supj] Ke,|P[;,][/:e,-)|}i/2. (6.14) 

|fc|>r 

Thus, by Fatou's Lemma and Cauchy-Schwarz inequality, 
E( hm sup||(/-P,)C/:P[^,,]e,-||) 

< liminfE({sup Ke'^I^Mf^>,)|}V2) 

k&Z'^\k\>r 

< liminf{E(sup |(e,|P[-,]t/:e,)|)}i/2. (6.15) 

\k\>r 

Now, by Fubini's Theorem, 

E(sup5^ Kefc|Pf-,j[/>,-)|) < E(j;supKefc|P[-,][/>,-)|) 

|fc|>r |fe|>r 

= 5]E(sup|(efc|P[-,][/:e,-)|) 

|fc|>r 

< CY e-"!'^--''!/^ < Ce"!-''!/^ ^ e-"l'=l/^ (6.16) 

|A:|>r |fc|>r 
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where the right hand side decays exponentially fast to zero as r — > oo. As a consequence, 
E( lim sup ||(/ - P,)C/:^P[^,,]e,-||) = 0, (6.17) 

SO that there exists a set C fl such that P(r2j) = 1 and for all uj G flj 

lim sup 11(1 -P,)t/:P[^_,]e,||=0. (6.18) 



In other words, P^^^ej G "^ppiUuj), Vcj G Then, (l = Cij^i^dilj is a set of probability 
one such that for all cj G Pj^^jCj G 'Kpp{Uuj), Vj G Z''. Hence, 



(6.19) 



Proof of Proposition 13.21 

By assumption, V = Ylk'^k^k satisfies Vfe = if \k\ > R, for some R > 0. Hence, by 
C auchy- S chwar z 



El 



|2p 



j \k\<R 



(6.20) 



i |fc|<_R 



since \{ej\U^P^(,-^ek)\ < 1- By the same steps as those performed in the previous proof, 
one gets that 

E(sup |||X|^[/:P[-,jV||) < oo if UirE(supKe,|C/:P[-,]e,)|) < oo. (6.21) 

j \k\<B. 



That the latter sum is finite follows from p.Sp . which ends the proof. 



7 One-Dimensional Localization 



In this section we prove Theorem [37 
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7.1 Basic Properties of the One-Dimensional Model 

The proof of Theorem 13.41 uses a number of tools, which are specific to the one-dimensional 
model, where the operators studied here can be considered as unitary analogs of Jacobi 
matrices. We start this section by briefiy presenting the properties which we will need. 
While we have a somewhat different point of view, much of this is related or equivalent to 
facts on CMV-matrices, which can be found throughout [491 150j . 

Due to the specific band structure of U^, the generalized eigenvectors can be studied 
using complex 2 x 2-transfer matrices. Specifically, generalized eigenvectors are solutions 
(not necessarily in Z^) of the eigenvalue equation 

U^il) = z^, il) = ^'ilJk(ik, (7.1) 



for z G C\{0} and characterized by the relations 



for all k £ Z,. Here the transfer matrices Tj : — > GL(2,C) are defined by 



t \ z 



+ ' 1 + e^(e-r,) 



(7.3) 



Note that det T^(6'^^, 6'2;,_^^) = e'-'^^'^^-^'^k+i) has modulus one and is independent of z. We 
have for any n G N 

t") =^^(^^("-i)'^^Vim)--^^(^o,^r) (^;-;) =T.(.,n) (^-) , (7.4) 

"^-'"-'^ -r.(r2„,r2„+i)-^---r,(r2,ri)-i f^-i") ^ r,(^,-n) f^:;") . (7.5) 



We also set T^(u;,0) = /. 

The transfer matrix formalism allows to introduce the Lyapunov exponent 7(2;); 

^ E(ln||T,(..n)||) 

n— »oo fi 

Positivity and continuity of the Lyapunov exponent for all values of z under the current 
assumptions was proven in [36]. A consequence of these properties of 7(2) is the following 
unitary version of Lemma 5.1 of |16j . 

Lemma 7.1. Assume that the random variables {Ok}keZ'^ satisfy i2. 8\) . then for each com- 
pact subset A of C there exist a = a(A) > and < 6 = (5(A) < 1 and C = C(A) < cxd 
such that 

E[||r^(w,n)7;||-^] < C7e-"" (7.7) 
for all z £ A, n > and unit vector v £ C'^. 
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We omit the proof of this lemma which is very simihar to the one given for the self- 
adjoint Anderson model in |16j, see Appendix A in [34J for details. 



We will frequently work with restrictions of the infinite matrices to discrete finite 
intervals [a, b] or half-lines [a, oo) and (—00, b], respectively. Here we slightly abuse notation 
and write, for example, [a, b] for [a, 6] HZ. To guarantee that the restrictions remain unitary 
we have to choose suitable boundary conditions. At each finite endpoint these boundary 
conditions can be labeled by a parameter in T. 

For a G Z and € T, the unitary operator S']f'°°'' on /^([a, cxd)) is constructed as follows: 
If a = 2n is even, let the 2x2 matrices Bi and B2 be defined as in ()2.2I) . then let [/]^"''°°^ 
be the unitary operator in P{[2n,oo)) found as the direct sum of identical i3i-blocks with 
blocks starting at 2n. On the other hand construct jj^^'"^^ starting with a single 1x1 



block e*'', then identical 52-blocks starting at 2n + 1. Now let s]q^'°°^ 
The operator slf"''°°^ on /^([2n,oo)) will have a band structure 



j-j[2n,oo)jj-[2n,oo) 
U e U o 



C'[2n,oo) 



re 



IT] 



V 



rt 
rt 



-rt 
-tr 



rt 
rt 



-rt 
-tr 



\ 



(7.8) 



The parameter r] can be thought of as a boundary condition at 2n. 

Similarly, Ue ig found as the direct sum of identical i?i-blocks with blocks 

starting at even indices, while Uo has identical ii?2-blocks starting at odd indices, 

with (C/i"°°'^"+^^)(2n + 1, 2n + 1) = e*'?. Thus, 



q(— oo,2n+l] 



rt 
rt 



-t^ 

-rt 
,2 

-rt 



rt 
rt 



-t' 

-rt 
,2 

-rt 



\ 



te'^i 
re^^ ) 



(7.9) 



To define °0:2n] slightly modify this construction, this time filling 

up and Ue with a 1 x 1-block e*'', respectively, yielding 



■,[2n+l,oo) 



Q(~oo,2n] 





rt 


-t" 




\ 






—rt 








rt 


r2 


rt 


-t' 




-t^ 


—rt 




-rt 


\ 






te''' 





(7.10) 
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while 



C'[2n+l,oo) 



re 
rt 



rp'2 

—rt 



rt 
rt 



\ 



-t' 
-rt 

-rt 



(7.11) 
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In 

a,b] 

Va,Vb 

have 



S, 



similar fashion, for integers — C)0<a<6<oowe can construct the unitary operator 
with rja boundary condition at a and 7]h boundary condition at 6, for example we 



c'[2n,2m] 



re 



rt 

rt r 

4-2 



t^ 
rt 

2 



rt r 



rt 

2 



v 



Finally, we define 

Tj[a,b] _ T^[a,b] o[a,b] 



-rt 



re 



M?2m 



(7.12) 



(7.13) 



where the diagonal operator d]^'^^ on /^([a,6]) is defined as in (j2.10p . Similarly, we define 



Uu}^r] ^ and U^^ri ' ^• 



and C/(i,..n°°'^' are characterized 



As before, the generalized eigenvectors of Uli^r^l^r^bi '^u),ri vuuj,-q 
by the relations (I7.2p are supplemented with appropriate relations to reflect the boundary 
conditions, see [34j for details. In the proof of Theorem 13.41 we only use ?? = 0, so for the 



rest of this section we write U, 



a,6] 



a;,0,0 



lj{-0O,b] fQj, 



simplicity, also frequently leaving the cj-dependence implicit. Other boundary conditions 
will be used later, see Section El 

In the following discussion we use the notation [/["'^l for general — oo < a < 6 < oo, i.e. 
we write U^-°°^°°\ C/["'~l and U^-°°'''^ for U, C/["'°°) and u'^-°°'''\ 

Another feature of the one-dimensional model is that Green's function G{k, I; z) can be 
expressed in terms of two generalized eigenfunctions to z which, separately at each endpoint 
a, 6, are square-summable or satisfy the boundary condition at the endpoint. 

For a solution 93 of (C/ — 2)99 = 0, we define ^ by 



</'2n 



rt 



ze 



rt \ I (P2n-1 



^2n 



A straightforward calculation shows that 93 is characterized by the relations 



^2k 
JP2k+l 



72k-li^2k 



'P2k-2 
^2k-l 



(7.14) 



(7.15) 
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for all A; G Z. Here the transfer matrices : — > GL(2,C) are defined by 

7;(0,r/)=ri(r/,0), 

where Tz{ri,6) is given by (j7.3p and T* denotes the transpose of T. 
For — oo < a < 6 < oo let 



(7.16) 



(7.17) 



To any z not in the spectrum of choose generalized eigenvectors (p"" and as follows: 
If a is finite, then ip"" is the unique solution of ([/["'°°) — z)yj" = with (p"'{a) = 1, i.e. 
a generalized eigenvector to z which satisfies the boundary condition at a. If a = — oo, 
then for ip"" we choose a non-trivial solution of (U — z)ip"' = 0, which is square-summable 
at — oo. In the latter case ip'^ is determined up to a constant (one can construct if from the 
tail of {U — z)~^eo at — oo and the fact that the transfer matrices (j7.3p have determinant 
of modulus one shows that there can't be two linearly independent solutions which are 
square-summable at — oo). Similar we choose ^p^ with prescribed boundary behavior at b. 

The following proposition gives an expression of the elements of G^""'^^ (z) in terms of ip"" 
and p^ and the corresponding ip"', ip^ defined as in (j7.14p . 



Proposition 7.1. For all finite k, I with a < k,l < b, if I = 2n or I = 2n + 1, 



G["'^](fc,/;z) 



cip\p'^ if k < I or k = I are even, 
cipftpl. if k > I or k = I are odd, 



(7.18) 



where ci 



Proof: A straightforward, if rather tedious, calculation shows that the matrix whose 
entries are given by the right hand side of ()7.18p is indeed the inverse of f/if '''^ — z, see |34] 
for details. I 

We conclude this section by proving the following lemma that is used later in the proof 
of Theorem [37 



Lemma 7.2. For z G C with < 
< G^(t, s) < oo such that 



2-E 



1 



2m) 



\t(p' 



2m-l 



+ (r 



ze 



1| < 1/2, a £ 
^<C,{t,s) 



2ml 



and s € (0, 1) there exists 
(7.19) 



f2m-l 
vim 



for all m > a + 2. 



Proof: First note that both '■p^m-i f2m independent of 92m and can not vanish 
simultaneously. This follows from the fact that the transfer matrices needed to construct 
them via (j7.2p from <p"'{a) = 1 only contain 6a, ■ ■ ■ , 02m-i- 

Therefore we have the following cases: 



29 



Case 1: (/?2m = 0, using that 

directly. 

Case 2: (^2m / 0- t^^is case 



^2m 



< 2|(/32rra-ll + '^\'P2m\^ bound folloWS 



27r 



d/i(6' 



2m J 



l^'^2m-l + (^-^e^^''")¥'2ml 



'2mh JO 



2tt 



(7.20) 



Let M = sup{|r — ze^^^mj . ^ [0, 27r],0 < ||z| - l| < 1/2} < oo and distinguish between 
two subcases; 



Case 2a: If t 



> 2M, it follows that 



t — \- (r — ze ^" 

f2m 



> 



V2m-1 



^2m 



On the other hand 



V2m-l\ * < (-o I J_^s| a is 



Using these estimates we conclude that 

/■27r 



< IJ_II2£ 2 +— ^ 



V^2m-l 
^2m 



(7.21) 



(7.22) 



(7.23) 



Case 2b: Assume that t 
such that for all /3 S C 

1 



I 

Jt 



/t 1(e±^^-/3)| 
see e.g. [4lj. Using this we get the bound 

f27r -j^ 



< 2M. For all < s < 1 there exists < C^(s) < cx) 

(7.24) 



< 



However, under the current assumption we have 



vim 



M 

< (4- + 2)^1^^, 



(7.25) 



(7.26) 



which gives the required bound. 
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7.2 Proof of Theorem [33] 



We now turn to the proof of Theorem 13.41 The proof is done in three steps: first we prove 
that it suffices to deal with even matrix elements of Green's function, i.e. that we only 
need to show (|3.15p for the case that k = 2n, I = 2m. This mainly serves the purpose of 
avoiding to cover four separate sub-cases. Then we show that proving the bound (j3.15p 
for an element (2n,2m) of {Uu; — z)~^ can be reduced to proving the same bound for the 
element (2n, 2m) of the resolvent of the finite volume operator f/P"-'^™-] ^t z. Finally, we 
prove that expectation of fractional moments of G^'^"'^™'' (2n, 2m; decays exponentially. 
Based on the Green's function formula (j7.18p this will be found by combining Lemma 17.21 
with Lemma l7. 11 



Step One: The following Lemma shows that the expectation of a fractional moment 
of any element G{k, I; z) can be reduced to the expectation of fractional moments of even 
matrix elements of G{z). This comes at the cost of enlarging the fractional exponent s due 
to the use of Holder's inequality. 

Lemma 7.3. Let s E (0,1/4) and k,l £ 7j such that \k — l\ > 4. Choose n,m £ Z such 
that k G {2n, 2n + 1}, / G {2m, 2m + 1}. There exists K{t, s, fj,) < oo such that 

1 

E[\G{k,l■z)n<K{t,s,^l) (E[|G(2r^ + 2i,2m-2j;z)|4^])'/^ (7.27) 

i,j=0 

for all z gC with < ||z| - 1| < 1/2. 

Proof: Using the definition of n,m and that |A; — /| > 4 one has |n — m| > 2. Since 
Lp^°° , defined by (j7.14p , satisfies (j7.15p and ()7.16p , a straightforward calculation shows that 

Using Theorem 17.11 along with (|7.15p it fohows that for k ^ {2m — 1, 2m} 

G(k,2m + l;z) = .^^ 1 [^2(1 + e-^(''2™-^2™-i) - ^ — —) - ze'^'^rn-i] 

rt(e"'2m-i -1/z) z 

G{k,2m;z) - t'^e''^^^r.-i-o^r.-2)G{k,2m - 2; z)} . (7.29) 

By Holder's inequality and (j7.24p it follows that for s G (0, 1/4) there exists < cl^\s, r) < 
oo such that 

E[|G(A:,2m + l;z)n 

< C7«(s, r) [mG{k, 2m- zt^\fl'' + (E[|G(A;, 2m - 2; z)^^])^/^) . (7.30) 
Similarly, using (|7.2p and (j7.3p we obtain 



= r(e^''2„+. _!/,){ l'^2n+2 + e -9^2™ }• (7-31) 
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Thus, for / ^ {2n,2n + 1}, s G (0, 1/2) and all z G C with < \\z\ - 1| < 1/2, 
E[|G(2n + l,/;z)|^] 

< Cf\s, r) ((E[|G(2n + 2, /; z)\^']f'^ + (E[|G(2n, I- z)\^-']f'^) . (7.32) 

It readily follows from ([7301) and ([7:321) that for |n - m| ^ {0, 1} and all s £ (0, 1/4) 

1 

E[|G(2n + l,2m + l;z)|"] < K';}\s,r) ^ (E[|G(2n + 2i, 2m - 2j; z)|^^])^/^ (7.33) 

i,j=0 

This proves the Lemma for the case A; = 2n + 1, / = 2m + 1. The other cases are more 
direct. H 

Step Two: Let [k,l] denote the interval [min{/c,/},max{A;,/}] . In what follows we 
show that the expectation of fractional moments of G{2n, 2m; z) can be reduced to that of 
^|[2n,2H|(2n,2m;z). 

Lemma 7.4. For s G (0, 1/3) and n,m £ 1^ with \n — m\ > 2, we have 

E[|G(2n,2m;z)n < C^{t, s){E[\G^^^'''^"'^\{2n,2m; z)^']^/'^ (7.34) 

for all z gC with < ||z| - 1| < 1/2. 

Proof: For definiteness, assume that m > n + 2, the case n > m + 2 being similar. 
Using the definition of [7<lf'^' ([7T3]) . we see that 



= ;7^-°o,2n-l] ^ ^[2n,oo) ^ ^ 

where is given by 

'{rt - t)e-*^2n-2^ 

_^2g — i6l2n-2 



(7.35) 



Kik,l) = { 



k = 2n-l,l = 2n 
k = 2n,l = 2n-l 
k = 2n, I = 2n 



(7.36) 



A; = 2n - 2, / = 2n - 1 
fc = 2n - 2, / = 2n 
(^2 _ ^)g-ie2„-i^ k = 2n-l,l = 2n-l 

(r^ - r)e"*^2n^ 

_i2g-i02n+i^ /c = 2n + 1,/ = 2n- 1 

(-rt + t)e-^^2n+i^ /c = 2n + l,/ = 2n 

0, otherwise. 

Denote G"(z) = G^-^'^""!] (z) G[2"'~)(z). By the first resolvent identity, we have 

Giz)-G-{z) = -G{z)KG-{z). (7.37) 
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Therefore, it follows for all m > n + 2 that 

G{2n,2m;z) = {I + t^e~''^^''-^G{2n,2n - 2; z) + rte~'^^^-^G{2n,2n - l; z) 

-(r^ - r)e''^^^G{2n, 2n; z) - {t - rt)e-''^^"+'G{2n, 2n + 1; z)} 
Gi^^'^){2n,2m;z). (7.38) 

A similar application of the first resolvent identity, this time to the difference G^'^"''°°\z) — 
Q[2n,2m] (^)^(;;[2m+i,oo) ^^l^^^ express G[2"'°°) (2n, 2m; z) in terms of G[2"'2m] (2n, 2m; z) 

and ultimately leads to 

G{2n,2m;z) = |l + t2e"*^2n-2(^(2n, 2n - 2; z) + rte"*^2n-i(^(2n, 2n - 1; z) 

-(r^ - r)e~*^2„(^(2n, 2n; z) - (t - rt)e"*^2n+i(^(2n, 2n + 1; z)} 
X |l - [(rt - t)G[2"'°°)(2m - 1, 2m; z) + 

(^2 _ r)G[2n,oo)(2m, 2m; z) + rtG[2"'°°)(2m + 1, 2m; z) 
_^2^[2n,oo) (2m + 2, 2m; z)] } G[2"-2m] (2„, 2m; z) . (7.39) 

If ^ and B denote the two {-j-factors on the right hand side of ()7.39p . then it follows from 
s < 1/2 and Theorem 13. II that 

E{\Af') < G, E{\Bf') < G (7.40) 

uniformly in |z| 7^ 1. Here we are using that Theorem 13.11 remains true with identical proof 
for the Green function of u!j^'°°\ An application of Holder's inequality to (|7.39p yields 

E[|G(2n,2m;z)n < (Ep|3^])i/=^(E[|5|3^])i/3 J^E[|G[2n,2m] (2^^ 2m; z)!^^]) , (7.41) 

which gives (I7.34p when combined with (I7.40p . I 

Step Three: Now that we have reduced the problem to dealing with fractional mo- 
ments of elements of the form G^^"'^*"" (2n, 2m), we show that expectations of such objects 
decay exponentially, in particular we show; 

Lemma 7.5. Assume that {9^}kei, are i.i.d. with probability measure d^{9) = T{9)d9, 
where r G L°°(T). There exist sq G (0, 1), < Gi < cxd, ai > such that 

E[|Gl[2"'2HI(2n,2m;z)|"o] < Gie~"il"-"l, (7.42) 

for all z with < ||z| — 1| < 1/2, and all m,n £ Z such that \m — n\ > 2. 

Proof: For m > n + 2, let ip"^^ and y^^™ be two solutions that satisfy the boundary 
conditions at 2n and 2m, respectively, such that ip2n = 1 and ip2m = 1- Using (j7.18p . we 
have 

pid2m 

G[2"'2™] (2n, 2m; z) = ~2m _ ~.n ^IZ- (7-43) 
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Since (/j^™" satisfies the boundary condition at 2m, it follows that 

~2m 



(7.44) 



Using this along with the definition of we obtain 
(^[2n,2m]^2n,2m;z) 



(7.45) 



Now, for s £ (0, 1) the expectation of the s-moment of Gt^"'^™'' (2n, 2m; z) is given by 



(7.46) 



E[|G[2"'2"^](2n,2m;z)n =i 



/ dll{B2m) 

Jo 



1 



^2n 



(7.47) 



where E is the expectation with respect to the random variables {d^} kez\{2m} ■ By Lemma 
17.21 we have 



E(|G[2"'2™l(2n,2m)n < C^(s,t)E 



Tju;,m-n) ( 

f2n 



(7.48) 



Using Lemma l7.H it follows that there exist ai > and sq G (0, 1) and cj^\so,t) < oo 
such that 



E[|G[2"'2H(2n,2m)|'^«] < cW(so, t)e-"i 



(7.49) 



for all z E C with < ||z| - 1| < 1/2, < e < 1/2 and m,n such that m - n > 2. 

In the case n > m + 2 let ip'^'^ and V'^"' be solutions of {U — z)ip = that satisfy the 
boundary condition at 2n and 2m, respectively, with — V'im = 1- Using (|7.18p . (j7.14p . 
we obtain that 



Q[2m,2n] (2n, 2m; z) 



(7.50) 



In order to bound the expectation we first integrate with respect to ^2m-i 
same procedure as before. H 



Proof of Theorem 13. 4t Without restriction we may assume |fc — /| > 4 (as (I3.15P for 

\k — l\ < 4 only requires the a-priori bound (|3.5|) ). Pick m,n ^ Z such that k E {2n, 2n+ 1}, 
I E {2m, 2m + 1} and |m — n| > 1. Thus using the results of Lemma 17.31 and Lemma 17.41 
there exists < K{t,s, fj,) < oo such that 



E[\G{k,l;z)\'] < K{t,s,^l) J2 (E[|G'[2«+2i.2m-2i]I(2„^2i,2m-2j;z; 

i,j=0 



|12snl/12 



(7.51) 
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Next the result of Lemma [731 gives that there exist soG(0,1),0<Ci<oo,q>0 
such that 

1 

E[\G{k,l;z)\'"] < C^^^K{t,so,n) e-"|2'^-2'"+2^+2il. (7.52) 

i,j=0 

Using the definition of m,n this yields (|3.15p . H 

8 Neumann Boundary Conditions 

In this section we study in detail the properties of one of the finite volume restrictions of 
the "free" unitary operator S introduced in Section TT-li This particular restricted operator 
will share many of the properties of selfadjoint Neumann Laplacians. We will therefore 
think of them as restrictions of S to finite volume with Neumann boundary conditions. 



8.1 Neumann Boundary Conditions for d 

Let e'l = r + it, L e N, and on ^^([o, 2L - 1]) define 



1 



-te''' 



i[0,2L-l] 



rt 
rt 



-t' 
-rt 



rt 



-rt r 



-rt 



rt 



r rt 
-rt r^ 
rt 



\ 



-t' 
—rt 

—rt re"' / 



which is a special case of the restrictions of Sq to finite intervals defined in Section [7.11 To 
characterize the spectrum of ^' define Afc = arccos(r2 — cos(/c7r/L)), /c = 0, 1, . . . , n, 

i.e. in particular A/, = and Aq = arccos(r2— t^), the latter coinciding with Aq as introduced 
in Section [2] and giving the band edges of as e^*'*'" . 

Proposition 8.1. The 2L eigenvalues of S^^''^^ are non degenerate and given by 



a{S 



= {e^^o,e*^^ = 1} U {e^'^^ : k 



1,..,L-1}. 



In particular, (t(S'|^'^^ C a{S). Moreover, 

ifo = {i,l,i,l,--- 
is an eigenvector associated with e^^° and 

^L = it,l,-i,-l,--- (-1)^+1)* 
is an eigenvector associated with 1. 



(8.3) 
(8.4) 



35 



To prove Proposition 18.11 we will use the transfer matrix formalism, i.e. solutions of 
Uijjip = zip are characterized by the relations 

t''^')=Ue^,,o-,^,)(^l'-'), (8.5) 

W2k+2 J V ^2fe / 

for all k £ 7j, where the transfer matrices : — > GL{2, C) are defined by (|7.3p . In the 
"free" case Sqi/j = e^'^ip, the transfer matrix takes the simple form 



T(A) := T,.(0,0) = ^ _ ^ _ J . (8,6) 

The following lemma will be used in the proof of Proposition 18.11 

Lemma 8.1. The vector (i, 1)* is an eigenvector ofT{\)^ if and only if X G {0, arccos(r^ — 
t2)}u{A : cos(A) =r^ -t^cos{kTT/L),k = 1,..,L- 1}. 

Proof: In order to simplify the analysis we distinguish between two cases: 

(i) (i, 1)* is an eigenvector of T(A). A straightforward calculation shows that this is 
only true for A S {0, arccos(r'^ — t^)} with corresponding eigenvalues {—1, 1} respectively. 

(ii) The second case is when (i, 1)* is an eigenvector of r(A)^ but not of r(A), i.e. 

i \ ( i 



T{\f ^ ; J = a ^ ; J , (8.7) 

while V = T{X){i, 1)* is linearly independent of (i, 1)*. Then it follows that 

r(A)S = r(A)^+i ( ^ ) = aT{X) l^i^=av. (8.8) 

Also since v and (i, 1)* are linearly independent it follows that r(A)^ = al. Finally using 
that the determinant of T(A) is one, we see that a} = \. 

Since the eigenvalues of T(A) are given by e='=*^''*^'^°'^((^' -cos(A))/t ) ^ have that a is an 
eigenvalue of r(A)^ if and only if e±*2Larccos{(r2-cos{A))t2) ^ -y^ -^ ^ 

A G {cos(A) = -t2cos(A;7r/L),A; = 1,..,2L - 1} 

= {cos(A) =r2-t2cos(/i:7r/L),fc = l,..,L}. (8.9) 

Combining the two cases gives the result. I 

Proof of Proposition I8.lt In light of the previous Lemma, it suffices to show that 
e*"^ G a(dff'^ ^^) if and only if (i, 1)* is an eigenvector of T{X)^. 

First it is not hard to see that e*'^ G a{S^^''^^ ^^) means the existence of "0 G ^^([0, 2L— 1]) 
such that 

( ^2^+1 ^ = r(A) ( ^^'"-i ^ , (8.10) 
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for m G {1, .., L — 2}, while 

Tpi \ ( -{r - e*(''"-^)) \ 



and 



J = ^^^-^ I _i(,_e^(.~A)) )• (8.12) 



Define V'-i sucli that 



Using (18. lip and that e^'^ = r + it, one can see that this definition is equivalent to having 
ip-i = iipQ. Similarly, defining such that 

\ rr^rxx ( 4'2L-3 



^p2L J V ^2L-2 y 

we see that ip2L = — ^V'2L-i- Also by definition, 

^'^>' ( t' ) ^ ( tr ) '"^> 

which shows the required assertion. The last two claims of Proposition 18.11 follow from the 
above together with r(0)(l,i)* = -(l,i)* and r(Ao)(l,i)* = (l,i)*. ■ 

We conclude this subsection with several remarks: 

(i) e*^''=™^('^'-*') G fT(5S'^-^"^^), while 6-*^"^^=°^ is not. 

(ii) For Ao = arccos(r^ — t^) = arccos(l — 2t^), sj^'^^ ^Vo = e'^^^^Jo has solutions with 
\Mk)\ = 1 for aU k G [0,2L - 1]. 

(iii) There is a gap between the upper edge of the spectrum of S'j^'^^ ^\ given by 
g«arccos(r -t )^ g^j^^ next closest eigenvalue. In particular, for any A; G {1, • • • , L — 1} we 
have 

|giarccos(r2-i2) _ arccos(r2-t2 cos{^fc/L)) | ^ ^2^^ _ cos{7rk/L)) = 2t^ siu^ {7:k /2L) 

using the property sin(x) > x{l — vr/4) if x G (0,7r/2). 

(iv) The previous remarks as well as the "Neumann-bracketing" property to be found in 
Section [TO] below will make the operators 5"]^'^^ a suitable tool when studying properties 
of finite volume restrictions of the unitary Anderson model near the upper edge e*^° of 
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the spectrum of Sq- To get the corresponding results also at the lower edge e of the 
spectrum of 5*0 one needs to modify the definition of S"]^'^^ by setting e^^ = r — it 
in (jS.ip . In this case we use, in particular, that the vector (— i, 1)* is an eigenvector of 
T(— Ao), leading to e~*'^° becoming an eigenvalue of the restricted operator. One gets 
properties similar to Proposition 18.11 and Remarks (i), (ii) and (iii) above. 



8.2 Neumann Boundary Conditions for d > 1 

For a box A := [2/i, 2mi - 1] x . . . x [21^, 2md - I] C Z'^ define 

= on ^'^^^l\[l,,2mj-l]) = l\A). (8.18) 

Note here that the discussion of Section 18.11 applies with obvious modifications to intervals 
of the form [2lj,2mj — 1] and integers Ij < nij. We will be particularly interested in the 
case of cubic boxes A^, := [—2L, 2L + I]"' for L G N. The spectrum of is given by the 
\Al\ = {4:L + 2)^ eigenvalues 

a(4^) = {n^tie''^^"^-.} (8.19) 

where 

kj E {0, 1, 2, • • • , 2L + 1}"^ for j = 1, . . . , d, 
(jj G {+1, -1} for fcj = 1, . . . , 2L, aj = 1 for kj £ {0, 2L + 1}. (8.20) 

Under the assumption ciarccos(r^ — t^) < vr, the upper edge of the spectrum of 5, 
gt(iarccos(r -t ) _ gidAo ^ belongs to cr{S^) and is non degenerate. An eigenvector corre- 
sponding to 6*'^'*'° is (fi^^ = iSiffo, whose components all have modulus one. 

Moreover, there exists cq, a numerical constant, such that 

dist(e^^^»,a(S;^-)\{e^^^n) > (8-21) 

Indeed, the closest eigenvalue to e''^^° is e'^^'^-^'^^o+^'^'^osir^-t^ cos{n/{2L+im ^ ^^^^^ -g 
degenerate, so that the distance (j8.2ip equals 

jgjdAo _ gj(((i-l)Ao+arccos{r2-t2 cos(7r/{2L+l)))) I _ jgiAo _ gi arccos(r2-t2 cos{7r/(2L+l))) | 

where cq = {tt{4 - 7r))2/8, see ([STTl) . 

For later study of spectral properties near the lower edge e"**^^" of the spectrum of S, 
we use the modified version of sj^'-"^"^^ from the fourth remark at the end of the previous 
subsection in the definition (j8.18p 
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9 The Feynman-Hellmann Formula 



The Feynmann-Hellmann formula provides, on the level of first order perturbation theory, 
the change of an isolated simple eigenvalue of a selfadjoint operator under an additive 
perturbation. Here we will need a corresponding result for multiplicative perturbations of 
unitary operators. We prove such a formula in an analytic framework, which will suffice 
for our purpose. 

Proposition 9.1. Let I cM be an open interval containing zero, "K he a separable Hilbert 
space and I B a U{a) an analytic map with values in the set of unitary operators on "K. 
Assume /3(0) £ S is an isolated simple eigenvalue of U{0) with normalized corresponding 
eigenvector (p{0) G "K. Then, there exists an open disc centered at of radius ao > 0, 
D{0,ao) C C, and two analytic maps D{0,ao) B a /3(a) S C and D{0,ao) B a ^ 
(/?(«) G !K such that 

U{a)ip{a) = (3{a)ip{a) Va G D{0,ao), 
dist{j3{a),(j{U{a) \ {/3(a)}) > 

||^(a)|| = 1 ifae D{0,ao)nI. (9.1) 
Moreover, for all a £ D{0,ao) H I , 

/?'(«) = {ip{a)\U'{aMa)). (9.2) 

Remark 9.1. i) For a given a, the last formula is of course true for any choice of nor- 
malized eigenvector ofU{a), corresponding to (5{a). 

a) If I B a ^ U{a) is analytic and takes its values in the set of unitary finite ma- 
trices, all its eigenvalues and spectral projectors admit analytic extensions in a complex 
neighborhood of I, even at the values of a where eigenvalues of U{a) may cross, see \4^ . 
Consequently, an analytic choice of normalized eigenvectors can be made for all a G I. 



Proof: By the general theory of analytic perturbations of operators, see e.g. [12], the 
operator U{a) admits an isolated simple eigenvalue /3(a), for small enough values of \a\, 
say in D{0, ao). Also, the analytic rank one spectral projector on /?(a), P{ct), given by the 
Riesz formula is analytic for all a £ D{0,ao). 
By definition, for all a £ D{0, uq), 

P{a)U{a) = U{a)P{a) = /3{a)P{a), (9.3) 

and since U{a) is unitary on the real axis, P{a) is self-adjoint for real a's. Now define the 
analytic operator M^(a) as the unique solution to the ODE 

W'{a) = [P'{a),P{a)]W{a), W{0) = I, a G D{0,ao). (9.4) 

It is a well known property ([42j) that the following intertwining property holds for all 
a G D(0, ao), 

P{a)W{a) = W{a)P{0). (9.5) 
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Note that l^(a) is unitary on the real axis, since its generator is easily seen to be anti 
self-adjoint there. We define an analytic vector by 

(p{a) = W{a)(p{0). (9.6) 

Identities (|9.5p and (|9.3p show that 

U{a)ip{a) = f3{a)ip{a) (9.7) 

and ip is normalized on the real axis, since W is unitary there. By diflFerentiation of the 
previous identity and application of P{a) to the result, we obtain 

P{a)U'{a)ip{a) + P{a)U{a)ip' (a) = 13' {a)P{a)^{a) + [5 {a) P {a) if' {a) (9.8) 

which reduces to 

P{a)U'{a)^{a) = I3'{a)^{a) (9.9) 
due to (j9.3p and (|9.5p . Since for all a G / H 15(0, ao) we can write 

P{a) = \^{a)){^{a)\ = VF(a)|<^(0))(<^(0)|H^-^(a), (9.10) 
the result follows. H 

As a specific application, let us consider the family of analytic unitary matrices 

U^{a) = D{a)S^ = diag{e-*"^^}5j^, (9.11) 

where is the Neumann restriction of S" to a d-dimensional box A introduced in Section [51 
a G M, and 0^ G T for all A; G A. U^{a) interpolates between and diagle"*^*}^^^, at 
a = and a = 1, respectively. Introducing the self-adjoint matrix = 'YlkeK^k \^k){^k\ 
on ^^(A), we can rewrite 

U^{a) = e-'^'^^S^, a G M. (9.12) 

Lemma 9.1. If e^'^^'^^ is a discrete non- degenerate eigenvalue of with normalized eigen- 
state f{0), then, for all a G M, there exist analytic eigenvalues e*'^*-"^ ofU^{a) with analytic 
normalized eigenvectors (p{a) such that 

Ae^AH ^ -ie^H»)y20j^\^ekMa))\\ (9.13) 

In particular, for all a G M, A'(a) = —J2keA^k \ {^k\v{oi))\'^ . 
Proof: 

The existence of analytic eigenvalues e*'*'*-"^ and analytic eigenvectors (p{a) of U^{a), 
a G M, follows from Proposition 19.11 and the remark following it. 
We compute 

(U^Yia) = -iHAU^ia), U^{0) = S^^ (9.14) 

and 

(vp(a)l -ii/^C/^(a)v9(a)) = -z ^ 0, |(efc|v.(a))|2e*^(°) (9.15) 

fcGA 



and apply (19. 2p . 
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10 Splitting Boxes by Neumann Boundary Con- 
ditions 

Throughout this section we will assume that boxes A C Z'^ are compatible with Neumann 
boundary conditions as defined in Section [U is given by (IS.lSh and 

= DS^ = diag(e-^^'=)5j^. (10.1) 

For notational simplicity we will assume in this section that the box A has a vertex at 
the origin, which does not cause a restriction. We first deal with dimension d = 1. 

Consider a one dimensional box Aq consisting of two disjoint adjacent boxes Ai and A2: 
Ai = [0,2^-1], A2 = [2/,2(Z + n)-l], Aq = [0, 2(/ + n) - 1], (10.2) 

with n, / > 2 (to avoid the special case 5{J'^'). We note that and U^^ © U^^ are both 
defined on /^([o, 2{l + n) - 1]). 

We want to show that the eigenvalues of U^'^ © U^^ are closer to the upper band 
edge e'^^o"*""^ of the almost sure spectrum S of [/ than those of . Recall here that in 
Section [3.2.31 we have assumed that \6k\ < a and Aq + a < vr. This is the analog of the well 
known property H^^ © H^^ < H^" , where is the Neumann restriction to a box A of 
the discrete Schrodinger operator. 

The following simple observation is the starting point of the analysis. Splitting a box 
by imposing Neumann boundary conditions is a rank one perturbation: 

Lemma 10.1. Let , j = 0, 1,2 be defined as above. Then, 

S^'^=S^^®S^' + \^){ip\ (10.3) 

where 

Tp = -te2i-2 - re2i-i - ire2i + ite2i+i 

f = t{-ie2i-i + 621) (10.4) 

The proof is an easy computation. This leads us to using the following fact about rank 
one perturbations of unitary operators which return a unitary operator. 

Lemma 10.2. Let U a unitary operator on a Hilbert space "K and f,g& 'K\ {0}. If 

V = U + \f){g\ (10.5) 
is unitary, then there exists (3 G (— 7r,7r] such that e*^ = 1 + {Ug\f) and 

V = e^W)(f\U, (10.6) 
where f = f/\\f\\. 
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Proof: 

The identity VV* = I implies that 

\U9){f\ + \f){Ug\ + \\gf\f){f\=0- (10-7) 
Applying this to / shows that Ug is proportional to /, so that 

U9 = (10.8) 



With this it follows from (fT03]) that 

V = {I + \f){Ug\)U 

= iI + {Ug,f)\f){f\)U 

= (/-|/)(/|+^|/)(/|)[/, (10.9) 

where = 1 + {Ug\f). Thus / — |/)(/| + = VU* is unitary, which shows that 

\fi\ = 1, i.e. fi = e'" for /? G (-^,7r], and / - |/)(/| + fi\f){f\ = e^'^l/X/l. ■ 

Taking into account the random phases, we apply the previous lemma to our case with 
U^o ^U^^ ®U^^ + \Di^){ip\ (10.10) 
and -0, ^ from (fTOll) . We compute = IIV'II = V^, i.e. ip = ip/\\ip\\ = and 

e^f^ = l + {{U^'®U^^)ip\Di^) = l + {{S^'®S^')ip\i;) = g-^ '^■-^^°«('-'-*') = e"^^" (10.11) 
so that 

f/Ao ^ ^-iXo\D4>){Di,\ jjAi ^ jjk2 ^ jj^-iXoMW ^ gA^ ^ (10.12) 

Let us introduce an analytic family of unitary operators defined in /, a complex neigh- 
borhood of [0, 1], by 

I3a^ U{a) = e-*°^ol^'^)WI U^^ C/^^ (10.13) 

such that [/(O) = U^^ © [/"^^ ^^d U{1) = U^'-\ By Lemma O we immediately get the 
following 

Proposition 10.1. Let a G I and e*"^*^") denote any analytic eigenvalue ofU{a), which is 
isolated except at a finite set of values of a. Then 

arg(A(l)) < arg(A(0)). (10.14) 

Remark 10.1. In other words, when a Neumann boundary condition is introduced to split 
Aq into Ai U A2, the eigenvalues ofU^i ® C^A2 '^'"g closer to e^^^°^^^ than those ofU\g. 
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Let us generalize now to dimension d > 1. Consider a box Aq of the form 

Ao = Ao(l) X A(2) X ...A(d) (10.15) 

where 

Ao(l) = [0,2(/ + n)-l], A(j) = [0,2/j-l], j = 2,--- ,d, (10.16) 

which we spht by a Neumann boundary condition perpendicular to the first axis as Aq = 
Ai U A2 with 

Ak = Afc(l) X A(2) X . . . A{d), k = 1,2 (10.17) 

and 

Ao(l) = Ai(l) U A2(l) = [0, 2/ - 1] U [21, 2{l + n) - 1]. (10.18) 
By the previous results, the corresponding operators U^'' , /c = 0, 1, 2 are related by 

[/Aq — g-iAo|Di/'>(DVi|(g)/®---(g)/ ^Ai _ (10.19) 

Applying Lemma |9. II again, with Hj^ replaced by the non negative operator \Dij)){Dip\ ® 
1® ■ ■ ■ ®1, shows that the spectra of and U^^ © U^'^ are related in the same way as 
in the one dimensional case, e.g. by (jl0.14p . Clearly, the splitting by Neumann boundary 
conditions can be done perpendicular to any of the d coordinate axes and can also be 
iterated. Thus we get the above form of spectral monotonicity also, for example, when 
spitting ?7^^ over the cube [— 2L, 2L + l]'^ into a direct sum of U{^^ for ((2L + '\-)/lY cubes 
Aj of sidelength 21. 



11 A Combes-Thomas Estimate 

Combes-Thomas bounds, originating from [19j, have become the standard tool in Schro- 
dinger operator theory to show exponential decay of eigenfunctions to eigenvalues which 
lie outside of the essential spectrum. They also provide a key step in localization proofs 
for random Schrodinger operators in the band edge regime, see e.g. [56] . Here we provide 
a Combes-Thomas type estimate for unitary operators with band structure. 

Let U be unitary on i^(Z'^). We say that U has band structure of width w > it can 
be written as 

U = D + with {ek\Dej) = 5kj{ek\Dek) and (efc|Oe-,) = if \j-k\>w. (11.1) 

Proposition 11.1 (Combes-Thomas type estimate). For a unitary operator U on l'^{7J^) 
with hand structure of width w, there exist < B < 00 which depends on U only, such that 
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Remark 11.1. i) The same result holds for U defined on a finite dimensional Hilbert space 
Z^(A), A C Z^, with constants independent of A. 

ii) Actually, our proof works more generally for hounded normal operators U with hand 
structure. Results of this type are known in the literature, e.g. f25^ . 



Proof: 

Let X = {xi, • ■ ■ ,Xd), where x„ is the self-adjoint multiphcation operator acting on e^, 
with k = {ki, . . . , kd), as x„efc = knCk and defined on its natural domain. We introduce the 
vector a = (ai, • • • ,ad) G M*^ and construct the self-adjoint operator 



e"^' acting as e^^e^ = e^'^Cfe (11.3) 



on 



= {V G l^Z") s.t. J2 KefclV')Pe"' < oo}. (11.4) 

Here ak = J^n'^nkn- Consider the operator 

:= e"^[/e-°^ = e^^Z^e""^ + e^^'Oe""^ = + 0„ (11.5) 
defined a priori on the dense set 

Co = {-0 G s.t. {ek\ip) = for \k\ large enough}. (11.6) 

The operator Ua is bounded because for any V' G cq 

Oa^P = Y. e'^(^-^^)(efc|Oe,-)(ej#)efc (11.7) 

\k—j\<w 

where the set F is finite and e"^^^-'^ < e'"'"'. Moreover, Dq. = D which shows that 
||C^a|| < Ci{a) < Ci < cxD on cq, for a in a bounded set. Similarly, 

\\U - UaW = \\0 - OaW < C2{a) < Cslal, (11.8) 

for \a\ small enough. From the resolvent identity, if z G p(U) and C3\a\/dist {z,a(U)) < 
1/2, 

{Ua - Z)-^ = {U- zy\l + {Ua -U){U- Z)-^)-^ (11.9) 

with 

\\{Ua-z)-^\\ < 2\\{U -z)-'^\\ < 2/dist(z,cj(C/)). (11.10) 
Finally, by the formula 

{Ua - z)-^ = e^'^U - z)-le-"^ (11.11) 
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we derive 



{j\e"^{U - zy'e-'^^k) = e-^'^'-^Hm - z)-'k) = {j\{U^ - z)-'k), (11.12) 
from which we get 

Choosing the components of a and their sign in such a way that IcKn] = |q:| > and 
d 

= ^oin{j-k)n > \a\\j - k\, (11.14) 

n=l 

we obtain the result, with |a| = dist (z, (t(?7))/4C3, and B = I/4C3. H 

12 The Genesis of Lifshits Tails 

After having introduced some tools in the previous two sections we will now start with 
the actual proof of Theorem 13.51 Throughout this proof we will focus on localization at 
the upper band edge e^i<^^o+°') of the almost sure spectrum S of [/^. The proof at the 
lower band edge is completely analogous. It uses the alternate form of Neumann boundary 
conditions discussed in Section [8] (setting e*'' = r — it va. (jS.ip rather than r + it) and a 
corresponding adjustment of the results on splitting boxes in Section [TOl 

We find it convenient to rotate the upper band edge of S to be identical with e*'*^" , the 
upper band edge of S. This is achieved by replacing the original by e~*"C/aj. In other 
words, setting Om = 2a we now assume 

supp/x C [0,^m] with G supp/i and 2d\Q + 9m < 2vr. (12.1) 

The latter means that S has the gap {e*'^ : dAo < ?? < 27r — (dAo + Om)}- 

As in earlier sections we will frequently drop the subscript uj from our notation. 

We will first establish a Lifshits tail estimate for the spectrum near the band edge e*'^'*'" . 
At the root of this is the following proposition which we prove by following the steps of 
Stollmann [56]. As in Section El for L e N we set A^ = [-2L, 2L + I]"'. 

Proposition 12.1. Let e^K^^^) ^ respectively e^'^^'^ , be the eigenvalue of largest argument of 
U^^ , respectively . Then X{U^^) < dXo and there exist b > and 7 > 0, independent 
of L and d, such that 

e^>^iU^^) - e"'^^\ < ^ <e-^^', (12.2) 

for L large enough. 
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Let us first give an easy Corollary of Lemma 19. li Recall that U^^ (a) is defined by 

dsn]). 

Lemma 12.1. Consider a fixed realization of U^'^{a) in dimension d > 1. Then the 
analytic continuation of any eigenvalue e*'*'^") ofU^^{a) is such that X{a) is non increasing. 
Consequently, je*'*'*^'^'' — e*'^(^^| is a non decreasing function of a > 0, as long as A(0) — 
X{a) < IT. 

Moreover, the eigenvalue e^'^^" of S^'^ and its analytic continuation e*'^°*^") satisfy 
Proof: 

The first statement follows from Lemma l9.ll and from 



1 2 

:2sin(A(a)-A(0))A'(a). (12.4) 



da 



The second statement makes use of the fact that the components of the eigenvector ip'^'^ 
all have equal modulus. H 

Recall the following standard large deviation estimate whose proof can be found, e.g., 
in Lemma 2.1.1 of 1561. 



Lemma 12.2. For non-trivial and non-negative i.i.d. random variables 9k and sq = 7o = 

-iln(E(e~^o)) > 0, we have 

^(^E^^^^oj <e-^°'^l. (12.5) 

Let us consider the small a behavior of e*^°^"^. 

Lemma 12.3. There exist ci > and c^ > 0, independent of d and L, such that for L 
sufficiently large 



da 



a=0 



< cia^L^, < a < ^. (12.6) 



Proof: 

Expanding e^'^^^i"^) in terms of q G M, we get that 

giAo(a) _ gidAo _ Q,_gjAo(0) ^ __r_e^^o(a) (12.7) 

da 2 d'^a 

for some < a < a. Next we use Cauchy's integral formula to bound 

1^' ""^ ^ = ^ / 7^^^^ = / —( (12.8) 

d^a 2m J\z-a\=r - ") 27ri J\-^^oL\=r \z - ay 
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for a small enough and suitable r > 0. Thus we need to control e*^"^^-* for z complex. Since 



U^^{a) - C/^^(0) = (e 



where 



/|| < gl^l^M _ 1^ 



we have, by the second resolvent identity, 

(U^Ha) - zr' = {S^^ - z)-\l - (e--^'^^ - " z)-'] 

for z a{U^Ha)) U (t(S^^). Hence, 

dist(z,cj(5^^)) > el°l^^^ - 1 =^ z£p{U^L{a)). 
Now (fOT]) says 



dist {e''^^\a{S^^)\{e''^^°}) > 5 



Thus, if \a\ < ao := . 



|A^|2/rf- 



I |z-e^°'^o| = 5/2} C p(C/^Ma))- 
We now take a G {—ao/2,ao/2) so that a < ao/2 and r = ao/2 so that 
{z I — a| = ao/2} C {z | jz] < a^} 



(12.9) 
(12.10) 
(12.11) 
(12.12) 

(12.13) 

(12.14) 
(12.15) 



and for such z's \e'^o{z) _ gidAoj < ^/2. Using 5/4 < ln(l + 5/2) < 5/2 if 5 < 1/2, one gets 
that if < a < 



a2 cp 



2 cPa 



P«Ao(a) 



< a' 



2 \aoJ 



We get the announced result with 



ci 



1286'lf 



C2 



(12.16) 



(12.17) 



t2co ' 32eM^ 
provided \Al\'^/'^ = (4L + 2)^ > 2*^00, i.e. for L sufficiently large. H 

Proof of Proposition \TTT\ Assume that \e'^(^^^) - e*^^o| < b/L^, with b to be de- 
termined later. Using the monotony in a (Lemma 112. ip and Lemma 112.31 we have for 
< a < C2/L'^ and L large enough 

d 



a 



da 



a=0 



< 



< 



< j^ + cia^L\ 



(12.18) 
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Dividing by a and then choosing a = Ci/L"^ such that C4 < C2 and C1C4 < so/2, we obtain 
that 



da 



a=0 



(12.19) 



< 6/C4 + So/2. 
Next we choose h such that 6/C4 < so/2 to get that 

< So- (12.20) 
Note that b is thus independent of d and L. On the other hand we have from (|12.3p that 

(12.21) 



da 



a=0 



da 



In probabihstic terms 



(4L + 2Y ^ 



(4L + 2)'' 



5^ ^fc < so > . 



(12.22) 



fceAi 



Finally an application of Lemma 112.21 ends the proof with 7 = 2'^7o. 



The Lifshits tail estimate of Proposition 1 1 2 . ll and the properties of the Neumann bound- 
ary conditions, Proposition 110.1] allow to prove the following result, which is based on an 
equivalent result for Schrodinger operators provided in [56] . 

Proposition 12.2. Let (3 G (0, 1). There exist finite positive constants 7, C and a sequence 
of positive integers Lj. with — > 00 such that for any k and any z £ C with 1 < |z| < 
2 and dXo — l/L^ < axgz < dXo, 



P(dist iz,aiU^^k)) < l/Lf) < Clf^-'^^^^e-^^^ 



d/3/2 



(12.23) 



Proof: Let (3 £ (0, 1) and b > the constant found in Proposition 112. ll Fix a constant 
C> 1. 

We claim that for each sufficiently large G N there exists Lfc E N which is a multiple 
of k and such that 



b 2 b 
To see this, note that ()12.24p is equivalent to 



Lk G 



i2k^/b)'/P,i2k^/b)'/^C'/^' 



(12.24) 



(12.25) 



As /? < 1, for k sufficiently large, this interval has length larger than k, allowing for a choice 
of Lfc as required. 
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We now show that (|12.24p holds for these L^. Split the box A^^ = [-2Lfc,2Lfc + 1]"^ 
into M = |ALj/|Afc| = {Lk/kf disjoint boxes as 

Kl, = ufLiAkU), (12.26) 

where A^ij) denotes the suitably translated box A^, + c{j). Consider now 

Un = U^'^'-^^ © U^'^'-'^^ © • • • © U^k(M) (12.27) 

on /^(Aij.) = ©^i^^(Afc(j)), where each [Z^^^-^) is provided with Neumann boundary con- 
ditions. Proposition 110.11 shows that passing from C/'^^fc to Un by introducing Neumann 
boundary conditions makes the eigenvalues come closer to the upper band edge, i.e. 

X{U^^k ) < A([/^'=(J")) for some jo S {1, 2, • • • , M}, (12.28) 

where e'^^^^^ denotes the eigenvalue of largest argument of U. As a consequence, taking 
the stochastic independence of the (U^''^^^) into account together with the relation ()12.24p 



< 2/L^j < P [^^^^K^^--^^-") _ e^'^Aoj < ^1^1 for g^rne jo 

< MP (^|e^A(c/(A.{i))) _ g.dAo| < 5/^2^ _ (;l2.29) 
By applying Proposition 112.11 to the box Afc(l) we see that the latter is bounded by 

^^e-i^' <CLf~^l^\--<"\ (12.30) 

Finally, since dist(z, cj(f/^-^ft )) < l/if for z such that Iz-e^'^^o] < l/zf implies |e*^(^'''^'= ) - 
gidAo| < we get the result. ■ 



13 Towards an Iterative Proof of Exponential De- 
cay 

The proof of Theorem [331 to be completed in Section [T^ will proceed as follows: To prove 
exponential decay of E(|G(/i;, I] z)Y^ we will join the two sites k and / by a chain of boxes of 
side length L. For a suitable choice of L and argz close to the edge of S, the Lifshits tail 
and Combes-Thomas estimates will show that the fractional moment of the finite volume 
Green function G^^\k,j;z) is small (think "less than one" even if this is only true up to 
some factors which can be controlled). Here G^^^ is the resolvent of a restriction of C/ to a 
box of side length L centered at k and j is a boundary site of this box. To turn this into a 
proof of exponential decay of the infinite volume Green function, we need two more tools: 
(i) a factorization of the infinite volume Green function into finite volume factors, often 
referred to as a geometric resolvent identity, (ii) a decoupling argument which allows to 
factorize the fractional moments in the geometric resolvent identity. These two remaining 
tools will be provided in this section. 

As explained at the beginning of Section [12] we will continue to focus on the localization 
proof at the upper band edge and continue to assume (jl2.ip . so that the upper edge of S 
is e^'^^o. 
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13.1 A Geometric Resolvent Identity 

Due to the specific structure of our operators (in particular their ergodicity with respect 
to translations by two) it is of advantage to cut up Z'^ into cubes of side length two. Thus, 
for n = (ni, . . . , n^) S let 

Cn ■■= [2ni, 2ni + 1] x . . . x [2nd, + 1] (13.1) 

and Xn ■= XCn characteristic function of C„. For L G N let 

Al= U Cn = [-2L,2L + lY. (13.2) 

\n\<L 

We will work with restrictions U^j'^ and U^j^ of U^j to Ki and its complement = 
Z"^ \ A/,. We choose C/^^ = D^^S^^, where S^'^ is the unitary Laplacian with Neumann 
boundary conditions from (|8.18p . In fact, the choice of boundary conditions is rather 
irrelevant as long as matrix elements are only affected near the boundary, e.g. we have 

Uf^Hj,k) = U4j,k) ifj,fcGAi„i. (13.3) 

Our definition of Neumann operators from Section [8] does not directly extend to operators 
on exterior domains such as A^, where the unitary Laplacian can not be defined as a tensor 
product of one-dimensional Laplacians. While it is possible to define Neumann boundary 
conditions directly for the d-dimensional operator, we choose a more simplistic approach 
and define 

U^l=pKu^P^i, (13.4) 

viewed as an operator on ^^(A^). Here P^l denotes the orthogonal projection onto ^^(A^). 

A'^ 

The price for our simplemindedness is that Ui^ ^ is not unitary. However, it is a contraction, 
i.e. ^ II < 1 and therefore a{Ucj{A'j^)) C {z £ C : \z\ < 1}, and it remains a band matrix 
whose entries satisfy, by definition, 

C/5(j,A;) = C/^(j,A;) ifj,keAl. (13.5) 

These properties will suffice for what we need in Section 

We will use what is often referred to as a geometric resolvent identity, relating the 
resolvents of ?7^, U^^ and Uuj^ . Following an argument which for the selfadjoint Anderson 
model is used in [5], we start by defining the boundary operator T^^^ through 

U^ = U^'- eU^^"^ +T^^\ (13.6) 

By the above construction of Uf^^ and Uuj'", in particular (113. 3p and (jl3.5p . the operator 
T^^^ has non- vanishing matrix-elements only near the boundary of A^, more specifically 

T^^^hx = XxT!,^^ = if|x|<L-lor|x|>L + 2 (13.7) 
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as well as 

XxT^''hy = ^ if|x-y|>2. (13.8) 

Also, the matrix-elements of T^,^'* are uniformly bounded in L and u. 

To keep the length of the following equations under control we will drop the arguments 
uj and z and for the rest of this section write 

G := {U^ - z)-^ (13.9) 

and 

G(^) := ([/Al e [/^^ - z)-l = (C/Al _ e ([/^L _ (13,10) 

We do a double-decoupling, once on and once on A^+i. Using the resolvent identity 
twice gives 

G = G^^) - G^^^T^^^G 

= G^L) _ ^ g^l)j.(l)(j^(l+i)q(l+i)^ (13,11) 

Observe that for y E Z"^ with \y\ > L+2 one has XoG^^\y = and xoG^^^T^^^G^^+^^Xy = 
0. Thus 

X^Gxy = xoG(^)r(^)Gr(^+i)G(^+i)x„ (13.12) 
which is the geometric resolvent identity to be used below. 

13.2 Decoupling of Fractional Moments 

The next result says that the fractional moment E(||xoGxj/||'^) can be decoupled along the 
boundary of A^. 

Proposition 13.1. For every s £ (0, 1/3) there exists a constant G = G{s) < oo such that 

iE(iixoGx,r) < c Yl mxoG^'-huin Yl mx^G^'-^'hyr) (13.13) 

u:\u\=L v':\v'\=L+2 

uniformly in z with 1 < |z| < 2, L G N and y G Z'^ with \y\ > L + 2. 

Proof of Proposition [TSTll From here on, the symbol C will denote a generic constant 
which may change from line to line but which depends on inessential quantities only. 
Define the boundary of Aj;, by 

OAl := {(x,y) eZ'^xZ'^iXxT^^^Xy /O} 

C {(x, y) e Z'^ X Z'^ : L < \x\ < L + 1, L < \y\ < L + l,\x - y\ < 1}. (13.14) 



51 



Expanding (jl3.12p over the boundaries of Kl and A^+i gives 

XoGxy = XoG(^)x«T(^^X.'Gx.T(^+')x.'G(^+'^Xy (13.15) 

(n,n') G Ml 
{v,v') £ dAL+i 

Taking fractional moments and also using that T^^^ and T*^^^^-* have uniformly bounded 
matrix-elements we get 

IE(||xoGx,ir) <C ^ (llxoG(^^X«inix.'Gx.||1|x.'G(^+'^X,ir) • (13.16) 

(n,n') G OAl 
{v,v') G 5Al+i 

Notice that the first and third factor in the expectation on the right are independent. 
Unfortunately, they are correlated via the middle factor. In order to decouple the factors 
we use a re-sampling argument, following a strategy developed in [3] and as a tool in 
the fractional moments approach to continuum Anderson-type models. For this, fix two 
pairs (n, n') G OAl and {v,v') G dAi+i- Let d := U C^' U C^, U C^,/ . In the resolvents G^^^ 
and we will re-sample the random variables On, n G 3- For this choose i.i.d. random 

variables {^n}ne3 with the same distribution as the 9n but independent from them. 

Noting that = X^^g^d e~*^"P„ (where Pn = {en,-)^n is the projection onto the 
canonical basis vector e^) we define the re-sampled -Dt^,cj '■= — D, where 

D := ^(e"*^" - e-^^")P„, (13.17) 
i.e. the variables {^n}n 

eg are replaced by the corresponding 6'„. Also define 

■■= D^,c.S^^ = ^^ - DS^^ (13.18) 
where ui^^ = U^^ , U^^^ = S^^ and 

G^^^:=(f^S-^)"'- (13-19) 
The resolvent identity yields 

q{L) ^ q{L) _ G^L)jjgKLQ{L) {U.2Q) 

and 

q{l+i) ^ ^(L+1) _ (13.21) 
We use this to bound the terms on the right of (113. 16p by 

lE(llxoG(^^X«inix«'Gx.inix.'G(^+'^Xyir) 

< EE [dlxoG^^^X^ir + llxoG^''^^4^G(^)x«ir)llx«'Gx.r 

(iix.'G(^+^)x.r + iix.'G(^+')^s^+^G(^+')x.r)" 

=: ^1 + ^2 + ^3 + ^4- (13.22) 
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Here we have argued that the above bound holds for arbitrary fixed values of the 9n- Thus 
it also holds after the average over these variables, denoted by E, is taken. Of the four terms 
Ai, . . . , A4, found by expanding the two sums in (jl3.22p . we will now find bounds for Ai, 
the one most easily handled, and A4, the most complicated one. Corresponding bounds 
for the two mixed terms A2 and ^3 can then be found by "interpolating" the provided 
arguments. 

Let E(. . . \3) denote the conditional expectation with respect to the tr-field generated 
by the family {0k}ki^3- Due to independence this means that 



E(X|3) 



X J] Tien)den. 



The re-sampled resolvents and G^^^^^ are independent of the variables 
by the rule 

E(X) = E(E(X|a)) 

for conditional expectations, 

Ai := EE[||xoG'(^^X«rilx«'G'x.inix.'G'(^+'^X,l 

= EE [\\xo&''^xurmxu'Gxv\m\\xv'&''+'^xyr 
< cEi\\xoG^^hu\nmxv'G^''^'hy\n- 



(13.23) 
^gg. Thus, 
(13.24) 



(13.25) 



In the last estimate we have used the bound E(||xu'Gxv||*|5) < C, e.g. Theorem 13.11 that 
the distribution of (^„)„gg is identical to the distribution of (6'„)„gg, and that XoG^^^Xu 
and Xv'G^^~^^^ Xy stochastically independent. This bound for Ai is of the form required 
in Proposition 113. ll 

We continue with A4, where we use ()13.24p again and then apply Holder's inequality 
to the conditional expectation: 



A, 



EE 



\Xo 



EE[E(|| ... 



^19)] 



< EE 



WxoG^'^^DS^^G^'^hufW^' 

xEiWxu'Gx.fV)'^' 
xE{\\x.'G^^^'^DS^^^^G^^+'hym)'/'' 



(13.26) 



We will now bound the three conditional expectations on the right separately. As s < 1/3, 
we can use Theorem 13.11 to bound the second factor by 



E{\\xu'Gxv\f'\dy^'<C<^ 
uniformly in {9k)k^2 ^^d z. 



(13.27) 
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To bound the first factor, we start from the definition of D to get 

||XoG'(^^^4^G(^^X«f'<C \\xoG^'^^Pit\\PiS^'G'^''^Xuf^ (13.28) 

^eanAi 

note here that xoG^^^Pe = for i A^. The unitary diagonal operator D^^ commutes with 
Pi and thus 

Xu\\ 

< 1 + |z|||P^G(^)x«II- (13.29) 

The re-sampled resolvent does not depend on the variables {0n)ne3- Thus, using 

Theorem 13.11 again, we get 

Ei\\xoG^^'^DS^'^G^'^^Xut\3) 

< c \\xo&''^Pif'na + \z\mG'^''hu\\f'\s) 

eeSDAL 

< C WxoG^^'^ Pef' ■ (13.30) 

In a similar way we find 

E(||x.'G(^+i)^5j,^+^G'(^+i)x,f 13) <G Y m4"^"G<^''^'^Xyf'. (13.31) 

By a calculation as in ([13:29]) we have HP^S^^+'G^^+^Xy II = ll^£Xy + zP^G'^^+'^'^XyW = 
IzlllP^G^^+^^Xs/ll- We conclude 

IE(llx.'G(^+'^i)5^^+^G(^+i)x,f 13) < C Y WPeG^'-^'hyf'- (13.32) 

Combining the bounds (|13.27p . (|13.30p and p3.32p into (|13.26p we arrive at 
< C Y ]EE(||xoG'(^)^'1l1l^'rG(^+'^Xyir) 

= c Y mxoG^'-^p.r) Y mPeG^'-^'hyr)- (13.33) 

Here it was also used that 3 has a fixed finite number of elements and thus (X^^gg x^*)^/^ < 

CY^j^gXj. The last identity uses that (0„) and (On) are identically distributed and that 

XoG^^^Pe and Pe>G^^+^'^Xy are stochastically independent. 

The bounds p3.25p . (|13.33p and related bounds for the mixed terms A2 and A3 combine 
via P3.22P and p3.16p to prove Proposition \Tm ■ 
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13.3 The Start of an Iteration 



We plan to use (|13.13p as the first step in an iterative argument, where the next step 
consists of applying (jl3.13p again, but this time with E(||xt;'G'(^^"-'^)xj/ll*) on the left hand 
side with v' as the new origin. However, before doing this we need to replace G*^^^^) with 
the original G, which can be done by reasoning similar to the decoupling argument of the 
previous section. 

Proposition 13.2. For every s £ (0, 1/3) there exists a constant C = C{s) < oo such that 

iE(iixoGx,r) < CL''-' nwxoG^'^huin E mx.'Gxy\n (13.34) 

\u\=L \x'\e{L+l,L+2} 

uniformly in z with 1 < |z| < 2, L G N and y G Z'^ with \y\ > L + 2. 

Proof: According to Theorem 113.11 we need a bound for EdlXi^'C-^^^^Xy II*) iii terms 
of fractional moments of the full Green function G for each fixed v' with ||f'||oo = L + 2. 
We start from the resolvent identity 

q{l+i) ^(l+i)^(l+i)^ ^-^3_35^ 

and expand 

^^^'''^= E X^T(^^'h^'. (13.36) 
(w,w')edAL+i 

Combining both yields 

iE(iix.'G(^+')x,ir) < mxv'Gxyin 

+ C Yl mXv'G^'^+'hn.rWXn.'GXyin. (13.37) 

(w,w')edAL+i 

With the goal of factorizing the expectation on the right we fix {w,w') £ dAi+i and re- 
sample over the variables 0„ for n G := Cy/ U Cw U Cw'- With independent random 
variables {On)^^^ independent from the On, but with identical distribution, define 



b 




(13.38) 










:= D^-b, 


(13.39) 




:= Du>^cjS = Uijj DS, 


(13.40) 


G 




(13.41) 



The resolvent identity 

G = G- GbSG (13.42) 
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yields 

< ]EE(||x.'G(^+^)x«,inix«,'Gxyr) 

+ tE{\\Xv'G^'^+^hu,r\\x^'GDSGxyr) 
=: B1 + B2, (13.43) 

where E denotes averaging over the variables 9n. Also writing E(. . . \3) for the conditional 
expectation with respect to the u-field generated by the variables {On)^^^ and arguing as 
in the previous section, one has 

B, = EE(Ei\\xv'G^''^'^Xn,\m\\Xn.'Gxy\ 

< GmXn^'Gxyin (13.44) 
Holder's inequality yields 

B2 < EE(E(||x.'G(^+^)x^fia)'/' 



X E{\\x..,GDSGxyt\d)'^') . (13.45) 
We have E(||xi;'G'^^+^)xw)|P*|3) < G and, by an argument as above, 

EiWxn^^GDSGxyfV) < CY,\\P,Gxyf'. (13.46) 

This leads to the bound 

B2 < GEE (Y^^WPeGxyl 

< G Yl IE(||x.'Gx,r). (13.47) 

Collecting (|13.43ll . p3.44|) and (|13.47l) into (I13.37|) . using that OAl+i has CL'^^^ elements, 
and ultimately applying Proposition 113.11 completes the proof of Proposition 113.21 I 



14 Proof of Band Edge Localization 

We finally have reached the point where all the main results of the previous three sections 
can be put together to prove Theorem 13.51 Specifically, we will use the Combes- Thomas- 
type bound of Proposition 111.11 the Lifshits tail estimate of Proposition 112.21 and the 
decoupling estimate in the form provided in Proposition 113.21 Also frequently enter will 
the a priori boundedness of fractional moments established in Theorem 13. 1[ 
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We will show the following fact which is equivalent to Theorem 13.51 (now in the nor- 
malization introduced at the beginning of Section [12] and only stated for the upper band 
edge): For < s < 1/3 there exist 5 > 0, a > and C < oo such that 

mXoG{z)xyr)<Ce-"\y\ (14.1) 

for all y G Z'^ and all z G C with 1/2 < \z\ < 2, \z\ / 1 and argz € [(IXq — d,dXo]. Here xo 
and Xy ai'e the characteristic functions of the cubes Cq and Cy introduced in Section 113.11 
Making the z-dependence explicit we write G{z) = (Uu) — z)~^ and G^^\z) = (f/^^ — z)~^ 
in this section. The bound (114. 1|) implies E(||xxG(2;)xy IH < Ce"^""^"^""" for arbitrary 
x,y ^l/ due to ergodicity. 

It suffices to prove (114. ip only for those z which in addition satisfy \z\ > 1. To see this, 
use the identity 

G*(z) = -UUU^ - z-^)/z = -U^G{l/z)/z, (14.2) 
which implies 

\\XxG{z)xy\\ = \\xyG*{z)x.\\ = \\xyUG{l/z)x,\\/\z\. (14.3) 
Inserting the partition J2y' Xy' and using that XyUXy' = for \y' — y| > 1 we conclude 

\\x.G{z)xy\\<^ \\Xy'G{l/z)x.\\. (14.4) 

\y'-y\<i 

This shows that (jl4.ip carries over to z with 1/2 < |z| < 1 once it has been proven for 
1 < l^l < 2, which is assumed for the remainder of this section. 

Proposition 112.2) shows that the probability that C/^^fe has an eigenvalue close to e*'^'*'" is 
small for the sequence found there. Combined with the Combes-Thomas bound Proposi- 
tion [TLT] this can be used to show smallness of the fractional moments E(||xoG'^^'=^(z)xu||*) 
on the right hand side of 1)13. 13p for values of z close to e*'^'*'°. 

Proposition 14.1. For any s G (0, 1/3) there exist a sequence of integers with — > oo, 
g > and C < oo such that 

E(||xoG(^'=)(^)x«r)<Ce-^''* (14.5) 

/or all k sufficiently large, any z G C such that 1 < |z| < 2 and argz G [dXo — Lj^ 2/{2+d) ^ dXo] 
and any n G Z'^ with \u\ = Lj^. 

Proof: 

Let 5l > 0, to be specified later. The Combes-Thomas estimate Proposition 111.11 states 
that there exists B > independent of L such that 

||xoG(^^(^)x«ll<fe-^'''^ (14.6) 
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for all z with dist(2;, o{U^^)) > 6l and all u with |n| = L. 

This takes care of the realizations uo such that the values of z are far enough from 
(t(C/^^). The Lifshits tail estimate takes care of the realizations where this is not the case, 
in the sense that such instances are very unlikely. 

We set 

= I dist (^,a(?7^^)) > 5l} and VLb = ^% = {uo\ dist (z,ct(C/^^)) < 5l}. (14.7) 
Making use of (jl4.6p . we can write by means of Holder's inequality 

mXoG^''\z)xur) = E(||xoG(^)(^)x«iri{.enG}) + IE(llxoG(^H^)x«iri{.en,}) 



2" 



+(E(||xoG'(^)(^)x«ir*)'/*(IE(l{.ef.,}))'/*', (14.8) 

with 1< t < 1/s and 1/t + l/t' = 1. Since E(l|^gj^}) = F{n) and E(||xoG'(^n^)x«ll''*) < C 
for < 1, by Theorem 13. H we get 

E{\\xoG'^'^Hz)xu\n < + C(P(dist (z,ct(C/^^)) < 6L)f/''. (14.9) 

To be useful for our purpose, this last quantity need to decay as L ^ oo, which requires 
L5l — > oo. On the other hand, we need 6l ^ for the probability that z is a distance 5l 
only away from a{U^^) to be very small, for suitable z in a neighborhood of the band edge 
gidAo^ In particular, this holds for the choice 6l = and any /? G (0, 1). 

More specifically, with this choice of 6l Proposition 112.21 yields the existence of a se- 
quence Lk with Lfc — > oo and positive 7 and C such that 

mXoG^'^Hz)xu\n < Ce-^^^^^Lf + CLf -^/^)/*'e-^^f (14.10) 

for all A;, 1 < \z\ < 2, dAo — ^/L^ < argz < cIXq and \u\ = L^. The choice /? = 2/(2 + d) 
leads to equal exponents of in the two exponentials on the right hand side of (|14.1Up . 
Choosing g < min(s-B, 7/^') and requiring k to be sufficiently large we can absorb the power 
terms in (jl4.10p into the exponentials and arrive at (jl4.5p H 

We proceed with the proof of (jl4.ip by fixing s £ (0, 1/3) and choosing the sequence 
Lk and g as in Proposition 114. 1[ We now also choose 6k = L^'^^^'^^'^\ 
Proposition 113.21 says that 

E{\\xoG{z)xy\n<cLi-' E mxoG^'^''\z)xu\n J2 mx.'Giz)xy\n 

\u\=Lk \x'\e{Lk+l,Lk+2} 

(14. 

if \y\ ^ Lk + 2. Let 1 < |z| < 2 with argz G [dXo — Sk, dXo]. This along with Proposition 114.1] 
imply, for k sufficiently large, 

E{\\xoGiz)xy\n < CLf-'^e-a^l^'"^"' E(||x.'G(z)x,ir). (14.12) 

\x'\e{Lk+l,Lk+2} 
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With the constant C from ()14.12p . fix L = for k sufficiently large such that 



b := CL2('^-i)e-^'^'^*'^"#{a;' e Z'^ : L + 1 < \x'\ < L + 2} < 1 (14.13) 
and get from (jl4.12p that 

E{\\xoG{z)xy\n<b rnax E{\\x,'G{z)xy\n. (14.14) 

Note that E{\\xx'G{z)xy\\^) = E{\\xoG{z)xy^x'\\''), which allows to iterate (|14.14|) . If 
x', x^'^^ , x^^^ , ... is one of the chains of sites obtained in this way, then the iteration may be 
continued as long as {x^^'^ —y\ > L + 2, i.e. at least — 1 times. For the last entry in the 
chain we use Theorem 13.11 to bound E(||Xa;0)G'(z)xj/||*) by C. In (jl4.14p this leads to the 
bound 

E{\\xoG{z)xy\n < (76^-1 = ^e^l^l. (14.15) 
Thus we have proven (jl4.ip with G = ^ and a = j"^^^ . 
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